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hip system prepared by 
Zediieses 2 Very structured preocadur>s for 
design-based comparisions and economic trade-off 
each and every stage cf the development of the desigr. The 
end result is a completes set of dasign requirements 
weeeneeO CONstruct the final preduct (this result sh 
meet the original design requirements!). 


This structured procedure has the form of the well-known 


Beecrguecpl sal" (Ref. 1], which is shown in Figure 1.1 for 
the propulsion machinery. This procedure #ncompasses the 
Bercolestop Gdesigr, including (but not limited to): 


Main Prepulsion System 
Shaft hersercwer 
Propeller RFM 
Specific fuel consumption and bunker capacity 


Space and weight objectives 


Bdevpeebility to ship configuration 


Steam-galley, dack, and heating systems 
Heating, ventilating, and air conditioning 
PrEcetohcang and ballastina 

Fresh water 
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Reno’ Shandiing 


SceC€ling engine and bridge taiametering certrol 
Garge handling gear 
Crane systens 

Flectzronic and Navigation Svstems 


Comminicaticn, ¢xtexrior and interior 
Radar 
Loran, Omega, etc., navigation aids 


Military electronics, sensors, command and control 
systems, weapons directors, tactical data systems, 


Price lecsrOPic Cou atermeasures. 


Mmueecnadce tc d=sign the vessel within economic constraints, 
pamcomusc =the life cycle cost as the function to cptimize, 
eprom @pcceei ve to define inivially the devices used <¢ 
provide the power for the vessel. They must provide nost for 
the least, in terms of the fuel that they expend and «he 
volume and weight that they possess. 

The selection of the main propulsion plant requires 
matching the power of the generating davice with a transmis- 
sion (in mest cases), the propulsor, other ship systems and 
weewnult.) &f One ligits the choices of propulsors to fixed- 
Feten EeEropellers and Controllable and rsversible—vcitch 
propellers, which are currently realistic for large vessels, 
the number of possible permutations remains quite sizable, 
2s shewn in Figure 1.2. The horsepower and RPM requirements 
cf mod¢rn warships require the use of a transmission to 
couple <he high RPM of an e¢concomical prime mover with the 


low BPM required at the propulsor for high prepulsive 
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efficiency. This transmission may be either mechanical or 
ezectrical. The mechanical reduction asar hes been 
exploited in the vast to minimize voluns/space requirements. 
When coupled to a reversible prime movar, such asa st¢an 
mee Oa 5 ae has keen proven to be quite acceptable 
Rowever, the electrical «ransmission has several attractive 


features, such as ease of propulsor speed and direction 
Gontrol, Preaxipi lity in the =he number and lLlocaticn of 
Con. ols, geeeecom of installation and machinery laycut, 
coupling of varicus prim? movers independently *o propulsor 
drives (as well as cther auxiliary applications) , and the 
ease of auxiliary thruster installation for maneuvering 
ecutrol. 
The choice of the propulsion. plant includes many impor- 

Gan factors, such as: 

oppo tiabi lity 

2. Maintainability 

3. Space and arrangement requirements 

4, Weight requirements 

Bomeyee Of fuel required 





6. Fuel consunmretion 

7. Fractional power and transient performance 
eee cere laticns with auxiliaries 

gs. Reversing capability 

10. Operating personnel 

17. Rating limitetions 
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Figure 1.2 Propulsion Plant Alternatives. 


These factors have gone through many cyclic changes 
based on the fluctuating economy. The recent emphasis on 
fuel ccests and availability has driven the fuel 
consumptionyprime mover efficiency to a much higher level of 
concern than it was previously given. The cost of acquisi- 
tion cf new vessels has skyrocketed, as well as the costs of 


maintenance and repair. These escalating costs are a chal- 
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Sepeauue-2cen tO its retirenent, elsO CcOlned me 2a crc b! 
C9styu 2S the driving fcsce that controls the desicn c= += 
mod2rn warship. 


the list of prime 


ce 
cx 
uo 


fe movers shown in Figure 1.2, 
eee cyere COS. Of the marines gas turbine is one of + 
a es The economy cf operation is maximized within a 
arrew range cf RPM's. This has been tolerated ir 
Citmee mp Oymen= of this device, since he cost of 
Meetazong an electric transmission to ae ee che GOal of 
ing the prime mover within its naximum economy range 
ema noe CEiset the lcss of efficiency Spherion aa raet by the gas 
turbine. This was due to the technology that existed within 
the field cf electrical contrels at th? time. The advent of 
puoetromac Eroqguency control of the output of the ganerator, 
Sedepernden= of RPM cr poles, and the ‘frequency centrol of 
the synchronous motor, independent of its supply frequency, 
makes the use of the electrical transmission much more 
attractive. Although the initial cost will be higher, the 
eis. Of OpStTabion, maintenance, ani repair of ‘the system 
Tatleeewialyeofiset the initial cost differential [ Ref. 2], 
ieee |, and ([ Ref. 4}. An analysis of the economic and 
engineering considerations involved with this mechanical vs. 
electrical drive selection was made by Eric Gott and S. 0O. 
evemccon | Ref. Sj in 1974; they clearly illustrated the 
Pacure cf the problem and possible solutions. ino heir 


paper, it was evident that the selection of some of the 


rossibly mere apprerriate forms of electric drives was 
discouraced due to the economics of the choices. High-power 
frequency-cecntrolling equipment was very expensive and 


bulky. Recent advances in high-power, solid-state equipment 
have now again premoted these alternatives into viabls 


contenticn. 
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A. GISTORICAL DEVELODMENT 


Powee S@ecil passenger launches, and tha Br 
experiments with small submersibles, to the present, the use 


@o Clec=rtic propulsion has n 


sv 
fu 
ww 


a slow evolution. The U.S. 
Wavy began with the Jupiter in 19171. It was a 5500 hp/shaf+ 
° 


Megeeteo-aie spec with wourd totor induction metors, which 


eum tec a 20-ycar lifespan (terminated by surface warfare 
ecteyaty ir 1943, while Serving as @n airsecraft carpier 


Vv 
(Langely)). This led to several major programs, summarized 


= AC drmve installed on five battl2 ships of 
30,000+hpd (1918-1922) 


* AC drivs applied to carriers Saratoga and Lexington of 
200,000hr (1920) 


PeoyaGazOnous ac -€rive irtroduced on four Coast Guard 


GUE ters (19°19) 


c+ 
an 
(D 
oe 
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mM 
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Qu 
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mw 
ct 


* Synchronous ac drive in 
160,000hp (1935) 


wesymer=renous ac drivé and full dynamic braking in 
cutters (1940) 

In the U.S.Navy Bureau of Ships Manual (BUSHIPS 
250-660=2), dated 1 March 1945, the text states that the 
Navy's inventory of ships that used electric drive (in addi- 
tion tc apnove) was "....about two hundred destroyer escorts 
and akou+ one hundred other ships including tankers, ca 
transperts, troop transports, and store ships." All of these 
applicaticns were with steam turbines (which were speed- 
adjustable and reversible) and eventually this design lest 
favor <tc the less expensive mechanical reduction ge 


systems which were somewhat lighter. The weight of these 


) 








umbe 
for spesd reduction caused the motors to 
£ 


BaGeeceevy7. The advent of the use of synchro He2Ors 

reduced the weight slightly since the motor experienced much 

fewer slip losses. AY g00d hastorical backgeound on the 
n 


poemimyeo-cO€S Of electric propulsion is found i 
Wenevensry Commander S.M. Robinsoa, USN [Ref. 6]. 

G@meesert applicaticns in the U.S. Navy are shown in 
apse I, 


TABLE I 
U.S. Navy Electric Drive Ships 


aye # Ships SHP # Mozors 
AS 2 15,000 1 

AS 2 qi 7520 8 
TARC 1 10,000 a 
TAGOS 1 1,600 Qu 


= also has 4-1200 hp thrusters 
fe 42lSo has 1-550 hp thruster 


BE. ELECTRIC PROPULSICS CONCEPTS 


Ameerea S£ current research in the field of rotating 


og 


electric equipment that shows considerable promise with 
regard to electric propulsion is the use of acyclic (homo- 
polar) superconductive DC devices. These devices produces 
excepticnally high tcrqu2 and horsepower in a very small 


volume and weight. This is of great value in advanced 


20 








UTZS70n EO (b= Censide req 
are dirs me (DC) 


devices. The ceategcry of AC devices is furthar subdivided 


nd Ses ae (AC) 


into inductive and synchronous devices. All three devices 
have keen used for marina prepulsion and each has its cwn 


inherert advantages and disadvantages. 


teueee Mechinery 


The DC electric drive system has been used 2xten- 
Sively in the past, primarily in low to medium power ranges 
(19000 - 6000 hpyshaft) with extremes of 400 hp/shaft and 
ie,o00 “Sho/shatt noted. This amploymert ef DC drives 
normally -utilizes 1-4 orim¢e movers powering one or more 
motors ver shaft (duplicating the system if nore than one 
Ebest exiSts) . 
Primary advantages of this system include: 
1. Mase of Control. Control being effected by varying 


the generator voltaga thiZzough field Vcontrol-=a process 


that lends itself well to remote location(s). 


Z- “ealtipie Control Stations. The ease of which the 
generator voltage and polarity of the motor conneéecticns 
(either field or_exciter) are varied allows multiple 
centrol stations to be designed for flexible ship 


Ger,c-O in tight maneuvering situations. 
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analysis). The 


Nas this dramatic functional 


Shewn in Figure 1.5. 
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Figure 1.4 Variation of DC Generator Weight 
with Rating. 


Spemacing targe is net realistic. Therefore, the remainder 


ef this paper will orly consider the AC systems. The 
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inherent losses associated with these DC systems though are 
essentially the same as those in the AC systems. Figures 
Vemeeibcugn 1.5 are from “Marine Engineering", Re diss 


Heme NQten, cds tor. 


Zemeee Machinery 


The cdevices ane buded jez the conventional AC 
machinery area fall into twe categories: induction and 
synchroncus. The principal difference betweer the two 
syszems is in the mctor epplication. The 2nduction motor 
system provides exceptional torque at low RPM and maximum 
slip (200%). This high torque response at low RPM makes the 
inducticn motor system preferable in some applications. For 
faVal abplications, this system is not as efficient as thse 


synchronous system by a few critical percentage points. Th? 
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Meeomeed De Woen the retor ané ths stator field. Thais siin 
Momeoeeeieca -C oho FPN end load, The remaining loecsss i= th 
mealiccren @Cro> arte the same as “he synchronous motor ir 
senetal. Another critical difference batweer <=he two 
systems is the power factor difference. ine Shey induc: 167 
MetOG, Sebe SO*tOor current ic induced from the stater field. 
Demecmeeretecsn 2isid “thet is created as a result cf this 


ent is the mechanism from which the mechanical rotatio 
erived. Since this Curren: is induced from the s ° 
D aVeuesent), it causes the current load seen by the 
supplying device to "lag" the supplying voltage ard there- 
fore cause the power factor (& measure of the amount of lag 

) +o be less than unity. A power factor of unity 
indicates that the load current and the load voltage are in 
phase and maximum real energy transfe> is occurring. The 
Beilaraon MOL5E is ncrmally operated with a lagging power 
BactO> @e epproximately 80% [ Ref. 6] and an overall effi- 
ciency cf approximately 90%. In the synchronous system, the 
Be@ge> fLaciOr is normally unity and is controllable. by the 
ercecrnal @o-Or rotor excitation which results in the synch- 
ronous system's power factor being adjustable and a unity 
power factor being easily obtained. Note here also that the 
rower factor in this system can be made to actually "lead" 
and thereby correct a system "lag" to other loads that might 
exist in the fully integrated system. The advantages cf AC 


systems in general are: 





ioh Efficiency. Losses in the synchronous system 
@esueaperox mately 6% (including excitation loss2s) as 
GCeémgraccad <co the induction motor System losses tha* 


2. Flexibility of Installation. The direct drive 


installad in tha vessel to minimize shaft 
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DE rot crs an a a 
reduction gearing and operated directiy co 


EZime mover. 


Power. The AC system allows 


Peers power from propglsion units tc be used for 


ae 
o 


wit reasonable upper li 
hp. The larger system's size is being reduced by 


taising the system voltage (which is not limit 


@ 
Qu 
os 
MG 
t 
> 
w 


commutation requirements of DC systems). The system 
Vermeages typically in use are from 2,300 to about 7,500 
volts. These vcltages are easily transformed for other 
M=2s teneaddition to prepulsion. 

The selection of the main propulsion system there- 
is not simple. The design environment must be specified 
he prime mover is selected. Onc2 these parameters are 


ed, the appropriate selection of transmission, control, 


Populser can be chosen. For the purpose of the propul- 


cf a vessel of approximately 8,000-ton displacement and 


ission of a nermal warship of the destroyer/cruiser 


iety using the marine gas turbine as the prime mover and 


Wete-=-cycle costs the desiaqn criteria, ene iselectien of the 


AC synchronous drive (possibly fully integrated) appears as 


@ webyearpropriate choics. Some of the design configura- 


<i9ns 
pertic 
(Ref. 
Feqgar 


moss ble within™ this choice, and the associated 
Imance of each, are shown in Figure 1.6 as presented in 
Wales Additional possible configurations are shown in 


ecmicimanc 8 . The end result of the implementation 


c+ these technologies would result in substantial savings 
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u 
feast ot] Tain propulsion in the life-cycle analysis. Suc} 
GeecoWPce2sOh if SHOwn in Figure 1.9, also from the WAYS: 
Bepo=t [ Ret. 7 j. 


Additional propulsion plant arra 
Semieneie Oo} and the total number possibl 
eoveyeoy) he designer's dimegination and creativity within 


Mmemcen- = 2ef s= the giuliforam 


SY eveLlVSeesstesye SS aoe 


ghee concep of I5Sses within ‘the notor and the 
generatcr ere comparative in magnitude and virtually iden- 
Peed a2 nature. The losses will, therefore, be discussed as 
a grcup which is applicable to both the motor and the gener- 
Bros The approach cf an energy balance is useful here to 
introduce these losses: 


{Energy into device] [Enerdy 1570 d-vico| 

( on the ie on the | = 

| mechanical side | | ¢lectrical side | 
| Increase of | | Ineresase of | ( Heat | 
mee croa Kinetics | { stored energy | |[generated| 
imecece-Jclamech. | + | in the elec.&mag.| + | within | 
| anergy in device | | €ields of device | | device | 


ale ale} 


The neat quantity on the right-hand side of the equa- 
tion is emall and positive; the majority of the eneray into 
the dewice is transmitted through the device as an output. 
As an example, when considering the motor, the energy into 
the device on the mechanical side is negative in sign and 
appreximately equal to the electric energy into the device 
mhackh 25 pesitive in sign. This amount must be balanced by 


the guantity on the right side of the equation. 
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* MeGnanical friction losses 
icss2s (windage) 

s and eddy-current losses 
* Resistance losses in the armature, and 
ses (in the case of synchronous machines). 
The ckiject for the electrical engineer is to minimize these 
goSse= te the greatest extent possible. Since it is impos- 
sible «+o eliminate all these losses, the problem that 
tremzins is the remeval of the resulting heat within <he 


device [ 


oO w 


@i. Si. The size and weight of the d¢vice is 
largely dependen+ on «he ability to remove this heat, since 
Piemicapecity >t the device is limited by the insulation usec 
in the device. The method of heat removal is not limited to 
a single moce, but is a combination of many modes of heat 
ereonstfer including: 

mecenduction: etator to casing 
EOLOr SOsstaccre (cieough eit: gap) 
EGror ©£o Shett 
shaf*= to casing 
shaft to external device (coupling, 


reduction gear, ¢tc.) 


weGenvec-10n (Doth forced and natural): 
casing to ambient 
etator toeint=tnal aia (Or gas) 
Teto= t6 antecrnal aie (or gas) 
° 


Shat= surface ¢ Sernal air (or gas) 





Shaft surface to external air 
aeRedzeticn: execanal surfaces $0 sinks 
imecstai sucsaces Eo casing 
These modes are rot all very controllable. tlserge 
aewnrts c& research have been devoted to mininizing these 
Puemm@al  2esietances. Efforts to improves convecticn hea 
Eradesoo have found that the restrictions on the rotor- 


Stator air gap have prevented exploitation cf this mcde as 
GmeevicSed as sarly as in 1926 by G.E. Luke (Ref. 9] and as 
ieeomes on 1979 by O.N. Kostikiv, et al. fret. 10]. The 
imprevement of the conduction heat transfer has been 
improved by minimizing the thermal resistance cf the insu- 
fe-2mg compcunds and the materials used in the construction 
cf the device, including installing thermally-conductive 
materials in the coil ends. 


The improvement of ccenvection heat transfer is still 


é¢ing examined by using liquids and gases to cool the rotor 
and the statcr of these devices. This vrocedure involves 
piping a liquid *hrough the stator to remove heat by forced 
convecticn, and utilizing a material with high thermal 
Copaucrwyaty directly edjacent to the piping to ccnduct it 
to the fiping where itis removed by the passing liquid. 
This me~hod is equally applicable to the rotor, provided the 
Pecowem c= tive rotating seals can bs overcome in order to 
Ghepmel the liquid from an external sink to the moving rotor 
and sbeck out. The advantage of doing this over standard 
cooling schemes is that the capacity of a device may be 
inezreased due to the lower internal temperatures that could 
be maintained (or the physical siza and weight of the device 


iss, 


could te reduced fer an existing capacity) as ¢hown ir 
Paguee 1.10 fom 2 40,000 hp, 180 RPM synchroncus motor. Th 


(wD 


use cf a liquid in a bouyancy driven closed-loop, called a 


themmosyrhen, is also possible. fhe amount of heat transfer 


a2 





Petqec-coOolingemethod in that it eli 

fOr, ECcatind seals. These séals can be weak points in the 
epee dQersom a relidbility stand point, especially if high 
ROMS are involved. The heat-pipe method typically involves a 
confined fluid acting as a two-phase medium fer heat 
jerensier, explceiting the phase change as a vehicle for 
substantial heat transfer in order to remove heat from th: 
BOPOreOr the Stator cr both. This method then transfers the 
eee wet-enes £5 the ambient air or to some other gas, liquid 
ee selsa heat exchanger. The configuration most suitable 


for use within a motcr/generator on the order of 40,000 hp 
would be either to ambient air orto a fresh water heat 
Sacmanecees <xXternal to the motor in the case of stator 
cooling, or +o internal forced air through extended surfaces 


Mn the case cf rotor cooling. 


C. PROBLEM DEFINSITICH 


Theyeurpose of this thesis is to analyze various methcds 
of cooling motors and generators of warship propulsion 
plants and to discuss the economic and reliability parspec- 
tives associated with each system. The next two chapters 
will discuss the ‘*theory and application of these two 
methods: liquid cooling of the rotor, both forced convection 


and the thermosyphon, and two-phase cooling of the rotor. 
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5 
Tevelopment Center, Annapolis, Maryland. Dt ses re cred 25 
and will be used for the analy 

sultsyreccmmendations herein may thereby be 
jiv evaluated, 

MO@mmraes DULPOSe OF this arbitrary coniiguraticn, ths 
model for analysis is as follows: 

1) The device will b2 a water-cooled frams, 

Syne =zorous AC qenerator, 25,118 KYA @ 3,600 RPM. 


2) A specific conduction bar shall be analyzed for a 
pemmecal load fer all cooling configurations. It has a 
length of 0.9144m and is 0.0116m squarewith 
a 4.763mm diameter hole bored through the lergth 
(kerein referred to as the tube length and tube 
diameter). The losses per bar in this configuration 
are epproximately 504. 
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Single-phase cooling has been extensively investicated 
mach regard to rotating references. Work was stimulated 
Seely ebweecrs Orinting industry, with the need to cooi 
Bei > ting 


lustry with the need +o cool high-temperature 
turbine Flades, which rotate at high RPM, and by the elec- 
eerveal “Beachine industry, where the need +o cool large, 
powerful generators krecame a necessity. Research in the 
study of flow and heat transfer in a rotating reference 
Syseemp “£5 andlicable to forced-convection flow, buoyant flow 
in thermcecsyrhons, and the single-phase flows of closed two- 
phase devices 

Mpemecometacal analysis of flow in a heated, horizental 
PuUDew wate reu- rotation by Morton [Ref. 11] was completed in 
oS 7. Although his analysis was limited to non-rotating 
flows, it gave considerable insight into buoyancy-induced 
secondary flows. His analysis was limited to laminar 
Gomvec-icn in horizcntal, uniformly-heated tubes at low 
Payleigh numbers (based on the temperature gradient along 
the pipe wall). The work was an extension of the study by 
Nusselt who ignored the gravitational (acceleratiocnal) 
erfects and, therefore, was independent of orientation. 
This study was a numerical solution to the governing equa- 
tions assuming small rates of heating and a notable varia- 


ceonmmorteedensity along the pipe length giving ris¢ to 


secondary flows. A Similar study was complet2d by Kays 
Rkecto V2 yin 1955. Morton studied this secondary flow 


driven by the bucyancy effects and the resultant temperature 
distribution by successive approximation in a powsr series 


(ee Is sami lar to the ~low ina pipe rotating about a 


5.8 





erc cular axis). Seceré 


D d 
BecOune: =©= the seccndary flom. 


The basic equations a 
Sepercomgen= 222 Bres@nted in Aprvendix A, and are the basic 
for many later studies. This development is adapted to «hs 
Beracing Case by the substitution of centripetal accaleza- 


S 
u 

M2on fOr the gravitational acceleration tern. 
e 


Maemee t= Gurarton of the cocling loop has =h>e¢ distinct 


Regions which must ke considered: the radial sections of 
tubing (oriented along the line of the centripetal acceler 

Ween), “che entrance regions of the axial sections, and the 
axial sections through the conductor bars (oriented parallel 
HOmeneeaxis of rotaticn). Additionally, it must be deter- 
ened whether the flow is lainar or turbulent The devel- 


pment will be reviewed chronologically within each topic 
With sample calculaticns made for pertinent correlations in 


Appendices Band cC as noted. 


A. HISTORICAL DEVELOPMENT OF ONCE-THROUGH COOLING 


ieeecne 1968 purlication, “Recent Advances in Hear 
fiomscer (ret, J3], Kreith made a survey of virtually all 
Premuor7. cc that date that had transpined inthe field of 
convection heat transfer in rotating systems. Has sece ion 
Gm concentmic cylinders and rcetating tubes summarized the 
works menticned herein, although principally evaluating the 
reports usirg gases as the fluid media. He presents the 
Sceeeptwor UELI2Z2ing a thermosyphon for the cocling of 
BOUO@s ,ecne »cheory cf which is antimately related +o the 
previcus work with single-phase heat transfer, which will be 
discuss2d later. Kreith also describes the three different 


sections (radial, entrance, axial) of the device which this 


study is addressing. These have been analyzed individually 
both expermentally ard theoretically. The inlet region has 
taememee! effect due +o the Coriolis force, which is a 
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Peeoneore ene ROS=Sty nisgber, or the swirl number (=: 
rocai of the Rossby number); the radial and axial tupe 
Meee Fevea, Wren a significant centrifu e 
€ me “NOtyceable hydrostatic effects ard, addi- 
tot 
ke 


s n 
ave, buoyant effects when a density dif 


ExXeezdang hic research to a vertical tube rotating 
about a perpendicular axis, Morris (Ref. 14] confirmed that 
the approximate solutions must includ@ second-order terms by 
studying low heating rates and using a method similar to 
Morton, but including tangential terms eas well as radial 
acceleration terms. He identified the three components of 
1) Forced convection due to the externally-generated 


pressure gradient, 


2) Gravitational buoyancy in the axial direction, and 
3) hOsee .oNa. buoyancy due =o centrifugal and Coriolis 
forces. 


Mae telative magnitude of the terms due to item 3) is char- 
acterized by the Rossby number. Ghosts ea ratio of the 
Mercia force to the Coriolis force. Although the Rossby 
number exists in the velocity and temperature fields for 
Eouetoneastorough the second order, it has no effect on the 
soluticns for heat transfer or flow resistance. Since at 
eneenehinal cotation rates and eccentricity that exist in 
BOtateng ebectric machinery, she centrifugal acceleration is 
much much greater than ‘the gravitational acceleration, the 
foraer wiil dominate. 

tue thsitr 19668 re 


Rodman Grr=ctlat Pares in the laminar region, Mori and 
+ 


A 


Ore vOnenecatetranste= inerctating, 


Nexeyema (Ref. 15] con ed their rasearch by evaluating 
° 


u 
steady, fully-develcped flow with analytical techniques. 
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Pe=noq. =he bcoundary-lavar concept, ney Gormelated “heist 
Ses ts with experimental werk by Trefethen. They qeve 
Peer etacacns ©O> Eee (Nu/sNuo} Kecia G5 & EWrceie2 <¢ 2 
angutar velocity, Reynolds, Prandtl, and rotational Eeyneldis 
numbers for fluids with Prandtl numbers above ard below 
Mite ye(t.c., GJ2S2s, iiquids, glycerol, and liguid metals). 
teeters Udy Coniirms the stabilizing effect of the secondery 
flow cn the cr 


m c 
W itical Reyrolds number and notes that laminar 
So8Utsons ere still valid when the Reynolds numb 
@rge (12,600 o §) =500 rad/s). 
eeond Leport on the subject of rotating racial 
eepes, Mori, ere i. (Ref. 16] evaluated secondary flow 
fects due to the Coriolis force by assuming a boundary 
yer aleng the wall. With experimental results and corre- 
tions from previous work, they were able to show that the 
et of secondary flow is less in the turbulent region 
ewlieianar region, although still on the order of a 
increase over the value without secordary flow. The 
ns they presented agree with their previcus rspor+ 
[Ref. 15] and with experimental results with air; they are 
presented as 2 ratio cf Nusselt numbers based on an exponen- 
Real Lunceicn f the Reynolds and Prandtl numbers and are 
given for awide range of Prandtl numbers including gases 
ecd liquids. 
Ito and Nanbu [Ref. 17] experimentally measured 
Eeectrcm factors £Or flow in rotating, straight pipes of 
Epewleaectoss~secticn, and compared their results +o those 
of Trefethen and Mori & Nakayama. Tiete empirical) correla= 
Pome mseoreoticticn factor ratios for two ranges of Kt (where 
Kt = J2yRe and J= Qdesv ) are shown in equations (2.1) 
Agee Za2yecOr the turbulent region, At = Ktcer. These are as 


follows: 


44 





O-942 Kee.282, 1¢ Kt <500, (229) 


th 
~N 
~ 
rh 
So 
il 


Oa g42 Ku0- 09, Ke 500. (22) 


c= ee co ose tach saectOr practically coincides 


a 
Bienes ctetsObery friction factot given by Blasius: 


Hot le@inar flow, Wwoach Ccecurs When Kt is botow Ktcr, the 


follewing ccrrelation applies: 


f/f 5 = 0.0883 K19-25(1 + 11.2 K1-9.325), (2.4) 
a 
where Kl = J¥Eeo. Equation (2.4) is velid EC= 
@ees02< Ki <t07 and (J/Re)<0.5. Cher st avert. Ons actor 
Hien GUT rotation is given by: 
Eo = 64/Re (22.5) 
a 
Medica rnem CCCUrS at? 
Recreate C7 J). 23% 103 (2.6) 
Cr 
Re =1.07 JO- 23*103 (227) 
er 
Poe cue Lange 28<¢ J <2*10->. Calculation of the pressure 


Gcopeeweuld be required in the detailed calculaticns for 
desuane Cf a cooling systen, especially of the closed, 


=FOtaeang-loop systen. 
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—— = ee > amen me ee 


4 
arc 
ot = 


wD 


mug, net 2... Pa2t ey eeataetyeced th 
BeezeCeotis-ics of water flowing through a partially-fi 


Bepeewi-ch was rotating about its own axis and 


Whee 
Ma2mce2necd at a constant heat flux. The also summarized 
the 


from Jravitational to centrifugal forcs is accempa- 

{2¢€S 25 previously discussed by Tevilor 
Oj. Miewalenors calculated @ czitical 
h the Nusselt number changes its func- 


ws how the 


Ref. 
ic te 
nal dependsnce on the Reynolds number. The paper aise 
heat transfer is accomplished within a layer 
Pp 


e 
around the reriphery of the pipe. A GorErelation for this 


=a: 
<= 


film thickness is provided, and experimental results confirm 
peo Valicity and illu 


strate how the increasing depth of «he 
Pinel ehe rotating Elpe, up to completely full, does not 
isfer bevond the critical depth value. 


enhance the heat tran 
ctational Reynolds numbers used in the 


The range cf the fr 
experimental work was between 4,000 and 20,000; the enhance- 
ment (Nu/Nu) was approximately up to 240%. 

Cliver [ Ref. 21] analyzed nacural convection in a 
MmeaeZzon-cl cube with no rotation, following <he work of 
Colburn { ket. 22}. Tis Stucy, walle poireting out several 
problems with the criginal research and basic criticisms 
pee coe Ove Marcanelli (Ref. 23], combines the work of many 
Saemore sto, Obtain a Correlation (equation (2.8) ) for the 
Nussel+ number in terms of the Graetz, Grashof, and Prandtl 
humbers and the L/D ratio. The paper also points out that 
the Grashof number must be evaluated using the temperature 
gradient alcng the pipe axis instead of the gradient from 
the pipe wall to the fluid mean temperatur>. 

ome = 1.75 (Gz 45.6% 10-** 
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(ee ee ee 23 (2584 


i 
The correlation is supposed to incorporate the #:£f:2 


c ° 
natural convection along with forced convection by vector 
° h 


mw 

Qu 
fu 
Jee 
c+ 


Bgnyect- he nctes that the velocity prori 


i 
Peed) 2468 45 critical in the analysis and that natural 
Bomecc=icn 3h the case of horizontal cubes skculd 
peeadent of L/D beyond the entry region. For this c 
eo=releticn is also ovided for “*he Nusselt num 


he L/D term: 


Nu eye = 1.75 {GZ 
ct 0.0083 Gr Pr 0.75 0,333, e 
( a cy } (227) 
The study of the Taylor vortices was continued by 
u 


Pattenden (Ref. 24] 


ebout its own axis. Due to problems in the slip-ring eppa- 


Sing a nOrLzontale (axial) tube rotating 


ratus used in the experimental work, the errors in the accu- 


E2aceme@cectrement of th temperature precluded specific 


c 


termination of exponents for the Reynolds number in the 
Compe lation for the heatetransfer-coefficient equaticn. 
alitative analvsis of the enhancament due to tube retation 
was mads, noting that the axial-flow effects were small 
GCommperec to the rotational effects. Taylor vortices were 
neither ccnfirmed, nor denied; the flow was sufficiently 
turbulent to mask the effects that the vortices may have 
had. Rotation rates were varied to a maximum of 4,000 RPM. 
Mori, e+ al. (Ret. 25) ~tnvestzcqated Eurther into 
the effects o£ buoyancy on forced-convective heat transfer 
in horizentall y-oriented, keated tub2s with laminar flow. 
His study extended the work of Graetz and Nusselt, and went 
Reyord =re limitaticn of small heat fluxes such as those 


SeGuemedmbyecolourn, McAdams, Martinellz, et al., Jackson, ¢t 


Gy 





made serge L/D ratio. Theit study involved a larg3 <smpera- 
ure difieréence between the wall and the fluid such as would 
Seecirewltn large heat flux. The change in viscosity with 
temperature is approximately covered with the (» /y. )o-14 
tee We Sicec= of natural convection 


a ReRa product; ces eae Oe = 
cable in the range of ReRa less than 3,000, but within this 
Tange the effect of natural convection does not exceed 
several vercent. The Rayleigh number is a function of the 
BeolletneOows> Of the diameter; therefore, the effect that 
matural convection makes in small-diameter charnels is 
small. This illustrates that in ths thin layer near the 
well, the temperature and velocit change abruptly and are 
ese uCurteree. =Tom that in the inner part of the fiuid. 
When the ReRa product is large, the secondary flow is strong 


eeGmeptnc veloc: .y distribution in the axial direction is 
e 


utterly different frem th2 Poiseuills flow assumed in the 
previous works. The Nusselt mber is then calculated by: 

Nu = = ZO 

u geq7 kts -t_) ( 
mieber che Scationary Nussel= number for constant heat~flux 


Q 
fu 
“” 
(D 
+> 
i) 


Nu 


Peal. (2214) 


The experiments were conducted with air and fluids with 
higher Prandtl numbers; they were cnly considering heating 
Emme Orevicactonel  ftield (not an acceleration field due to 
EOtac 1h)’. teens case, the layer near the wall becomes 


thinner and the Nusselt number increases with the increasing 
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i) 
tid 
OQ 
oO 
if 
i} 
a 
‘ 


ReRa. The local Nusselt number for the experinenx: 


eeeen (2e= ait) is given by: 


Nu = 0.61 (ReRa)%- 2 (141.87 (ReRa) 0-2) (2.12) 
This study was extended into the turbulent region, where the 
Regmeest font of the temperature distribution was found to 
Pteereene che Girecticn of gravity slightiy. LPech2sS cass, 
the effect cf the secondary flow was not as significant as 
in the case oz laminar flow. The Nusselt numbers calculated 
within this regime were in good agreement with these of 


eorpuca [Ret. 22], where: 
Nu = 0.0204 Re9-8, (25743) 


Mori and Nakayama used the theoretical analysis of 

their previous research [Ref. 26] and [Ref. 27] on straigh* 
PeesesOrcaesng about a parallel axis [Ref. 28}. Analyzing 

e body forces driving the secondary flow caused by density 
feereceneces 20 the centrifugal field and the Coriolis force 
with regard to the flicw, they used fundamental principles to 
characterize the flow and temperature fields. This was done 
assuming an effective secondary flow due to buoyancy and 
using the rotational Reynolds number. This would indicate a 
Beped csWergqence frem the results of Morris [Ref. 14] and 
eowe that Coriolis effects cannot be ignored. The paper 


Gave correction facters as a function of Reynolds, reta- 


tional keynolds, and rotational Rayleigh numbers, and 
concludsd that «he Coriclis effects diminish with increasing 
eccentricity (higher g-fields) to a negligible alue when 


the ReRa product is greater than 10,000,000. 
In 1968, Nakayama { Bef. 29] further analyzed a hori- 
Peele Stte Gur pipe rotating about a parallel axis. He 


assumed, as previcusly [Ref. 15], Siew an s er ect ives 
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Secemgcoy LlLOW in tuliy-developed con 
a 


fu 
}4- 
c 
| 
Oo 
Fe | 
M 
Oo 
th th 
th 
}- 
Oo 
2! 
rey] 
os 


erature fields existed, and included body 
Bee ore Correlations he obtained for both = 
i 


c 
emectric=10n factors and the ratio of Nuss2ii 


Wet d ove= 2a l2raqe range of Prandtl nunbers for both liquids 

and gases, as well as for awide range of Reynclds and 

Grashcf numbers. The correlation for the Nusselt number for 
d eu 


(140.014 /7(Rey p2-5) 1/6)} (2712) 


where Beeera<ce/is (Grrpro.6)—12/13, theca term is the 
PeesorOte sne=stia force to body force. Inertia fore] is 
represented by Re22/13, and the remai 
represent the body force. The body f 
S@orgwebhas fo0rse in the case of the inlet re 
metvons, Of £ke tedial a@zms, or the centrifuga 
@ase Ge the axial section. The numerical result of his 
ectrelations for tke model presented in Chapter I is 
meeluded ~nm Appendix B for beth turbulert and laminar 
irauids. 

Steqwasth, e¢: al. (Ref. 30] further evaluated the 


effect that bouyancy-induced secondary flow would have on 


t- 


aminar flow in a heated, horizontal pipe without rotation. 
is work shewed clearly that secondary flows did exis* due 
to density variations and that the heat transfer was 
equaeced by themas a functicn of the Grashof and Prandtl 
numbers. This effect is somewhat negated in *he presence of 
ext=eme acceleration as that in the periphery of a rotor, 
but remains nonetheless, provided the turbulence is below 


the Critical Reynolds value. 
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Peer eco eossine 1 Ots as tO the Ststctuse sf the flc 
end temperatuce fields to simplity the governing equ nS, 
Boome 22a Pore: Ss oresented a numerical eslution «+o tha cass 
Gemeene SOta-ing cylindrical tube [Ref. 31] with laminar, 
Bore y-aqevelcped tlow. Data were compared to the thecretical 
Eesuits for the case of air, water, and glycerol. The 
eo= 2 J2220n obtained is in terms oe che 


Piyeeeten-neynoids-OSranct] preduct as is given by: 


Nu = Nu, (0.262 (RaRePr) 9°173}, (2.15) 


The value for the test model is calculated and presented ir 
Appendix B. 

tepherson again studied this parallel, TO=e =210- 
Fipe, heat-transfer problem for fully-developed turbulent 
£VON He compared his experimental results and correlation 
meehethe Carlier work of Morris and Woods [Ref. 32} and they 
compared favorably within the turbulent entrance region. He 
Herecs ha: he rotational bucyancy was net a strong factor 
in the secondary~flow in the fully-developed region. 
Because of this, His results, when compared against 
Nakayama's results assuming a strong ¢ffect [Ref. 29], did 


not compare well. His correlation is list2d below: 


Nu = 0.0071 Re0-88 J0-.023 (2o10) 
and is comvared with previous correlations for the test 
model in Appendix B. 

3. Emtrance Regions 
In 1969, the Institution of Mechanical Engineers, 


Thermodynamics and Fluid Mechanics Group, sponsored 42 sympo- 
siun, MPeeticesibject Of heat “2>anslezr and fluid flow in 


electrical machines. Included in the presentations was a 
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pee oy 6 Deviacscn (ef. 33], which described «4 fuli-siz3 
os 


Beeeha-Or=rotor <es Eiqd52¢6= Pevealuating the heat=<t2rancisr 
Peecectsrisetics cf hydrogen cooling. His testing snowed 
thaz enhancement in the turbulent flow regimes was net as 
Mmagh as theorstically predicted. It was hypothesized that 
he teason Was the inability cf the hydraulic-diameter 


Gemeet. ©C sCoupe ely account for the flow conditions in «he 
mon=—circular ducts. The tests were able to give a good 
turbulent correlation with theoretical prediction when ths 
cooling scheme was mcdified by shortening the axial path. 
This effect, although not explained, could have been due to 
the Coriolis swirl effect in the inlet region enhancing the 
heat transfer above that to be obtained in the fully- 
develcped region 

The effects of the entry length, especially with low 
Prandtl numbers such as those of gases, are very difficult 
+9 eliminate and the Coriolis:'effects on the heat transfer 
may be noted even at high Reynolds numbers. As the 
Feynolds-Rayleigh product increases above about 1,000,000, 
there is a tendency for the amount of enhancement to 
dig@inish. This is attributed to the turbulent effects over- 
iding the secondary flow enhancement. The effect of rota- 
ticn on heat transfer is to eénhance it in the entrance 
regions ty the secondary flow due to the Coriolis effects 


ana in ere €ully-develcped, horizontal-tube regions by buoy- 


ancy effects. A sinilar survey of the technology was made 
by Petukhov and Polyakov in the USSR with much the seme 
results [Ref. 34]. The heat-transtfer problem in the 


eme=anees region of atube, where the Coriolis terms doni- 
nate, is addressed by Morris and Woods (Ref. 32]. 
Correlaticns are presented for air and are applicable for 
cther gases as well. These correlations, fOr both <{he 
laminar and the turbulent cases, are presented as functions 


GE the proaduet cof the Reynolds number and the 
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Cc 
ie@errUsueirsh WOTK iS ERequiTed in the cass G2 liqnics, 
Fi ard "akevane 


cS Nicene Cle Ho 


y of the stata-ocf-the-art tect 


also presented a surve 
ETegarding the heat-transfer characteristics of 
mepes ena ducts [{Ref. 35}. The questions of th 
errect and its length were emphasized as beinag generally 
unanswered, but they gave a general guideline for deter- 
Meneng =he e€xcent of the sntrance région as approximately 20 
*imes the diameter of the pipe. They also verified correla- 


tions cf orevious works for helium and water. 


4. Combinations cf 


sntcaq in Ehe 1969 Institution of Mechanical 


es) ne) 
4 
oO 
Yn 


osium was a study by Lambrecht (Ref. 36], 
he problems associated with water cooling 
S paper summarized the theoretical considera- 
as reviewed his previous work and the work of 
esearchers (Neidhoeffer and Ingenieure). He 
perior ccoling properties of water and the size 
eed Weight reduction, along with the improved efficiency of 
water-ccocled machines. He listed values for the optimun 
cooling duct size based on the heat-transfer and 
electrical-loss characteristics for various gasés. The 
method used for optimum duct calculation was based on the 
principle that the pressure drop in the duct decreases with 
increasing duct size and electrical losses increase with 
@eczeasing conductor area (increasing duct siz¢). A suit- 
@ble Calculation of this type is necessary in any final 
Gesign fcr cooling of these devices by any method devised 
SOrweTerecEecne Single- or the two-vhase approach. 
In 1970, Sakamoto and Fukui measured heat- and 


Si s5-troererer cocftficients for air and oil, specificelly for 
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Boome neccoling rorating slectric machines [Ref. 37]. Noting 
meaeetn] MOS: effective and Economic method to increase the 
PPP Ge oor o> Slectsical sachines was *c Lnorcvs “he 


~~ 8 Sm ~ wy ree er ale oj 


‘ 


EeOiing C= tho insulating material (sspecially in th> rotor 
eonductors end drums), a gecmetry similar to a typical ro 

meas tabricated and suitable coolants were tested. The 
Geuwecctive heat transfer was characterized as a functicn of 
with the H/d parameter held 


Seomb/a, ke, srr, Pr variables, 
mene Shape factor Gf eccentricity H/d is 


eons =an=. when ac 

sufficiently large, the centrifugal acceleration may be 
assumed to be a unifeorm acceleration ( 2H) acting en the 
Pioewax2s, Which 2s similar to the force gravity exerts on a 
mo==ZeON-al tube. With a lew Reynolds number, the empirical 


Cos cedation was in the form: 


Nu = Nu, (1+#C (Rar) ° 1 76z) ™ (2. 17) 


where the values of the constants were: 


© ="0-03 tee 0.75 Nn = 473 
and where Nu = 1.86(RePrd/L) °° 33 for large L/d values, 
and Nu = 2.67 Gz9-33 for small L/d values. 


This was valid within +10% to -17% for Reynolds numbers 
retween 162 and 2,700 and rotational Rayleigh numbers to 
20,000,000. It was also found that the length of the 
entrance region seemed to have no affect on the correlation, 
which weuld imply that a secondary flow existed in both 
eonfmgureticns. An Inceresting poin® ih this paper is that 
+he rotational Grashcf number (hence the rotational Rayleigh 
number) was calculated using the temperature difference 
bretween the wall and the fluid instead of the axial tempera- 
ture gradient. 

woods and Mcrris investigated laminar flow in the 
rotor windirgs of directly-cooled, electrical machines in 

e 


(eevee tekees Sejeeand noted that if water was used as the 
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Seomeme =n 4320 Of hydrogen (which is currently widely 
[eee enh Str2eC=sncy improvement of >0.54 could be octained. 
Suyen the. the 0.5% cf a large, say 1,000 MW, generat 
the only consideration made (and neglecting th 
weight reduction 


é Z 
the extended operating life ro 
@us to the cooler 0 y 


and 
operation), they proceeded + 
meinmdamentels o£ the froblen. They surveyed th 
ag 


Merks Cirectly app 


1) Morris {Ref. 14] 


iat asrecr~, low reotaticn ra 


series-expansion technique 


a 
2 and low heating rates), and 


2) Mori and Nakayama { Ref. 28] assumed a secondary flow 
(which was claimed to be valid fer high rotation rates), and 


used en integral-type analysis. 


woods and Méorris claimed both were inadequate due to the 
Booe- te J@nhc 1mposed Ey the nature of their solutions. They 
attempted a versatile and "exact" sclution for laminar flow 
in the fully-deaveloped regicn by solving the governing equa- 
tions with a numerical procedure. Numerical solutions wer? 
Beesented fOr both the friction factors and the Nusselt 
numbers. Their analytical values wer2 compared with #xperi- 
mental results and to correlations of previous works. 
Discrepancies were explained and the difficulty in obtaining 
fully-develcp2ed conditions in the experiments of their work 
and previous works was emphasized. The axial density- 
variaticn effects were noted, also, as a potential for error 
Pfeecheis enalysis at high values of retational Rayleigh 
tumbers. 

Nakéyana and Fuzicka evaluated the generator 
problem, where reasonably large radii of otors (approxi- 
Ba-ely)1 @) and high RPM (3600) creates ae centrifugal 
eccelerationr affects [Ref. 39]. This acceleration causes 
TioMmecenemicugal buUCyancy term to be Significant. It also 


See 





eteeeeces wo Cozlolis 22m in the radial tube seqnerts end 
-_ = = + « bs - == ~ a —_ ~~ —— ~ 
Sueved anls > recions. SiS oer oneen co=relations fos tas 
usssit funhecs, 


i-Tac  Oceraee oser fol Gach 5 


v t i 
B-ceons (fadial, curved inlet, and axial), and compared the 
ef 1 s 


BesuitsS tO ¢xpetimental data. Their correlations are listed 
Eslows 
mate = 262 RO-9-33 (radial) (2.18) 
Beytie = 1.5 RO~o.3 ices) 21nd) (219) 
EF/fo = 14.0 RO-9- 8 (axial) (2. 20) 
Nu = Re0-8 5-0-4 F0.014 (RE sRo2-5) 0-124} (radial) (2.20) 


Mein r=cert confirms the results of Nakayama's correlation 
omme@e- 14) ) .tor the axial séctions. 

The application of the technology vor ro 
cooling schemes is net limited to the printirg industry and 
electric machine area; extensive research has be d 
eomrnis) field of study by the gas turbine industry. The 
Advisory Group for Aerospace Research and Development 
(AGARLD) held a meeting in September, 1977, at Ankara, 
Turkey, to discuss High Temperature Problems in Gas Turbine 
Engines. Numerous papers were presented on the subject 
ancluding a paper by W.D. Morris en flow and heat transfer 
in rotating coolant channels [{fef. 40]. He used a selection 


of experimental results to illustrate the influence of 





Oo 
ct 
ions are modified to inclu 
SeCOrsOlusehocces aid 
r constant properties, these 2q 
M eae Gistences sufficiently downstream fr: 
J where the axial gradients of the velocity 
Gible, the Gliminaticn of th pressure-gradient 

2 tne 


a radial and tangential momentun s 
Coriclis terms to vanish identically as a source for the 


equations caus? 


creation ct secondary cross-stream flow; see also (Ref. 14}. 


ingehe entry region, the Coriolis terms interact 


x 
rw) 
ct 
on 


the developing axial velocity and create seccndary 


Berpendtcular tc the axial direction, even with the 


ch 
O 
Oo 
x 
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ch 
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wee cepe-ty ccndition. Finally, 2m the § EulLly= 
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Ge Of the Sotaticon manifests 
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rm) 
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ped region, the 
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nly as a cross-stream pressure distribution if buoy- 


|4- 
nate 

Qa nN <s 
wo 


ry} 
ry 

ct 
kG hk 


2tms ar2 neglected. If the density is allowed to vary 


x 


. 
3 


temperature, the centripetal acceleration terms need 
Siige@ese included as the Coriolis terms do not affect +his 
acticn. This cross-stream, buoyancy-induced flow gives rise 
te greatly-enhanced keat transfer and an attendant increase 
in flow resistance. The SrtEsces Of Ene eccentricity, with 
megardscc the Coriolis terms, hes @iecle effect, provided 
Biassthe ratio of the radius (H) to the diameter of the tube 
(d) is greater than 5 (only very small radius of rotation 
will have any noticeatle effects). 

Marto (Ref. 41] lists the papers presented at the 
14th Symposium of the International Centre for Heat and Mass 
Meas=ft== (TCHAT) held in Dubrovnik, Yugoslavia, in 1982. 
Papers of particular interest to this study included a paper 


by Jchnscn and Morris [Ref. 4&2], wherein “he concept of 
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ang itS influencs was shown 70 
Bee cerolily avydrtesta tic. Phesc© ork also confirmed the 
SG@ecct | CL COL101li1s acceleration with regard te e¢corndary 
=lows, especially in sections with relatively small L/d 
ratios. No quantitative recommendations were vresanted. 


MGrTi PReft. 25] states, when considering turbulen<« 


{> 


eloe, that it is not necessary to consider the influence of 
Eaoyancy on hsat transfer. The secondary flow was found to 
suppress the turbulence level when the turbulence at the 
inlet region was high and an empirical formula for the crit- 
ical Reynolds number, in terms of tha Reynolds number and 


Rayleigh number product, was calculated: 


Re = 128(ReRa)o- 25, (25-22) 
Cr 
Contrary *o this, when the turbulence at the inlet region 
was low, the critical Reynolds number was higher (7,700 vs. 
2,000), and heating decreased the critcal value. The net 
resul+ was that when the ReRa value was high, the secondary 
flow caused by buoyancy makes the critical value of Reynolds 
number tend toward the same value, whether the turbulence in 
tee inlet is Ligh or low. In the case of low turbulence in 
the inlet, the value for the critical Reynolds number was 


given by: 


Re = Re /(1+0.14 ReRa*10-5). (2.23) 
Cr cr, 

Mie Graph Sshewn as Figure 62 in the Krieth article [Ref. 13] 
is shewn below as Figure 2.1. This graph shows the heat- 
Beanee=t as 2 functicn of both the RPM and the flow-through 


Reynolds number. It is easy to see the flow's transitional 
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The transiticn to turbulent flow is suppressed by 
rot2ztion and has been qualitatively noted by numerous 
authors. An exact correlation for transition to turbulent 
flow in a retating-reference frame with heating is net known 
Z Ome ars % From the available literature, it is clear that 


the transition follews a path from the laminar correla*ions 
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EG ene turbulent correlations. Pia S Pech, on) the Cice. of 
Becmeeme-- 1 TAaference, has a “dip as shewn quelite+ively 
gy 2 pict of the Nusselt number against the Reynolds numben, 
Paeqgure 2.2, Whereas in the non-rotating case, no “dip” 
eccurs. 
eo A 
es 
f ¢ 
fo 
y Za 
(iT 
fe 
» 
\ where; 
ke 2S & 
Re, 
S 
Z 
-Turbulent -Laminar Region - 
Region - 


(All Re, > 2300) 





REM 


Re, > Res > Re. Reg 


F2gure 2.2 Deans; ttonel Chie racrceristics of 
Rotating Systems. 


6. Summary 


Eased upon the inform nh reviewed in the 1litera- 


atio 
ture, it may be concluded that the three regions within the 
motor/generator application (shown in Figure 1.11) hav? very 


G@itterent heat-transfer characteristics. 


Wem tne radial section cf the coolant path, the Coriolis 
EQuces cend tO domanate the heat transfer and friction. 
This effect is to enhance the heat transfer and increase the 
pressure drcep in these sections. The heat-transfter enhance- 
ment is substantial and could be effactively exploited for 
she heet tejection from the device by axial fcrced-air 
Gonmecticn Of che substantially-hollow, synchronous rotor. 
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resul+ is enhanced heat-transfer and frrere 
ee 


effects, but to a lesser degr 


Byrn the “@anding bars, the centripetal forces doin 
mese £Ohces result in hydrostatic effects which 
a 


Y 
Siew would pe affected by the Coriolis forces a 


He] 
ai | 
it 
: (p - 
(n 
yp C 
t- 
(t 
4 


b 

some degree of seconcary flew due to buoyancy. The 
cf secondary flow is questionable and is somewhat affe 
by the degree of turbulence and the amount enhancemen 
less specific. 

The enhancement in the laminar-flow heat transfer 
has been reported by numerous authors. This can be seén in 
Figure Da. Bio tgea tha Nakayama Rei. 294 send the 
moods-Norris {[Ref. 31] correlations for Nusselt number fall 

m © valte of 48711 and the @ushbu- 


lent value cf the Dittus-Boeslt 


between the classical la 
er correlation. They havé the 
Same trend throughout the range of RPM shown and differ by 
only a few percentage voints. Avuparently, the effect of 
rotation in laminar flow enhances the heat transfer because 
cf intense secondary flows du2 to centripetal ferces. An 
ghceresting point is that, at the higher RPM values, the 
corrélaticns of Nakayama, Woods-Morris, Stephenson, and 
Dittus-Bcetler all are exceptionally close together. 
Mrewhneat—transfLer predictions for turbulent, ftully- 
develcped, parallel flows are shown in Figure 2.4. This 
fagure illustrates that che resultant heat transfer is some- 
mere 2m the vicinity of the classical turbulent correlation 


cf Dittus-Bcelter for the current model. 
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The heat transfer is ¢ither scmewhat enhanced (on the crdar 
Ww 4 


s 
Sess in the lo Zee SGlera va iat ROt 29 jeje 35 
Ss 


ei a “ Sy SS ee ter as , = S = = Peat einen ne Sale eta 
Se@ reat Suppressed ty the high centripetal acceleratis 
£ ~ = } “re 4 <= o Gat 
Pomectie as Sho¥n by the cotrelaticn of Stephenson { Ref. 43}. 


t of the heat transfer in the Nakayama 
J in the lower RPM range is assumed *o 
u 


be due to the resi PCO ctetase orf eces etroi ce entrance 
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Be HISTORICAL DEVELOPMENT OF THERMOSYPHON COOLING 


The thermosychon is a device tha*+ creates the motion of 
Bae wO2King fluid sclely by buoyancy forces. Tt has a 


regicn where heat i 


“" 


Supple to =hes working fluid sand 4 


region where heat i 


Nn 


rejected. The less buoyant, ccoler 
Posewocne ct =h2> fluid moves in the direction of che accelera- 
meenmvectou, While the opposite is occurring for the warmer 
BoEtaOn. Tha +hermcesyphon also has athird region, the 
Coupes ag ©eqioa, With problems of fluid shear, mixing and 
entrainment. Fiaure 2.5 illustrates the conceptual arrange- 
ments of the three basic thermosyphon systems. 


Kreith (Ref. 13] mentioned in his review that the ther- 


mosyrhon tye. , a liguid-=-filled container) had a great 
Peeeretal £o: cooling of slecttic rotatin machinery. In 
motors and generators, the use of a thermosypvhon would 


preclud> the necessity of using rotating seals, which have 
weegqem@etesceliabilgjty in applications of high RPM and volu- 
Beieric flow. In 1973, Japikse [Ref. 44] gave a thorough 
review of the literature regarding thermosyphon technology. 
Beginning with a reference +o the Davies and Morris paper 
eeet. 95), which will be discussed in detail later, Japikse 


categorizes +hermosyfhons as follows: 


61 





1} tha Sera of hot nG2 ey =e = ar 3 ~ 1a wa ~— AN 
we od ac =~ — = ws OMOle! ees {e2.< aS CA rmcie Se eee he = 
Ress <clow), 
- k= = 2 = moh a ” =~ - 7 - = 7 vy os 
ee eC me 6S Reet threansts= {purely netusel cony-ic- 


) 
P2on cre@ixed natural and forced convection), 


3) the number and types of phases present, and 


4) the nature of the body force present (qravitaticonal 
aesOra.20nal). 
Gress troad estegqcries are the oq On gates 


m 
4 
wD 
< 6 
ot 
rH 
[ie 


C 
prescribed circulating fluid system driven by thermal tuoy- 
eacy forces." 

The use of the closed thermesyphon in rotating refer- 
ences, suchas in gas turbine cooling, has been studied in 
detail. In 1965, Bayley and Lock [Ref. 46] axperimentally 
investigated the perfermance cf closed thermosypvhons. They 
uzilized the theoretical develooment of the open thermosy- 
phon, modified the theory, and experimentally verified their 
Lesuilts. itewes Shown that the critical operating p 


ters include: ; 


ijaeeLenaqrca—to-D2emeter Ratio. This retio controls the 
characteristics of the flow and determines the extent of the 


eoupling region. 


2) Heated-Length-to-Cooled-Length Ratio. As this ratio 
memdeesoward zero, the analysis is identical to the oren 
thermcsyrphon correlations. TaLS Teseticts =he use of thse 
analysis for the clcsed thermosyphon to Lh/Lce ratio to be 
not much qreater than one. 


Bympccuplang Regacn. This is the region that compli- 
cates the analysis and causes the Prandtl number to affect 
the heat +t 


eonduceson, convecticn, and mixing. These are regarded as 


ransfer. Rte a<1sts sn theres distinct modes; 


G2 





feos ata eCela: Coupling ususlly consists of 2 conmbinatior: 


s so 
ing, heat-exchanger fins, home 
n 


ic cool-down Bees 
er ° 


steam *ubes 5g } ovens, internal combustion en 
eooland, and envirecnmental control in the space vehicls 
Sogaam, aS “Weil as many others aentioned in the Japikse 


Im 1971, Bayley and Maetin f{Ref. 47] also rev 
Searcoee-cene=abe technoiogy, with particular emphasis on 
gas-turbine applicaticns. They studied both the cpen and 
closed =hermosyphon systems. The us¢ of an open thermosy- 
phon system in a rotating reference, such as in gas turbine 
S2orimd, as well as in the cooling of electric, rotating- 
machinery area, has problems regarding the fluid selection 
and its saturation temperature and pressure. In the open, 
rotating reference system, the region that receives the heat 
is under high acceleration forces with high resultant pres- 
sures. Rsee-ehe hacural £low due co buoyancy occurs, the 
heated fluid moves frem a high pressure to a relatively low 
pressur2 area where a phase¢ change tO vapor may occur, 
possikly even explosively, blocking the flow. They intro- 
duced the idea of using two-phase cooling within the closed 
thermosyvyphon as away of exploiting this phenomenon, by 
using the high heat fluxes associated with the latent heats 
of evaporation and condensaticn and the lower temperature 
gradients associated with the phase change processes. This 
also has an added benefit cf reduced weight over the 
Single-phase liquid system. They noted that further 
Eesearen in the area of the Coriolis force effects on the 
heat-*ransfer chacteristics in a rotating reference needed 


to be accomplished. 
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Byene 2he tech 


Moen d 22 =ecnnology CE the tharmosyphon to electric 
BeCa™nes,, fcrris and Davies [Ref. 45] adapted the principles 
ame c | c.csed =het=mcsyphon to = cle Pec CO hee = 


Gould be placed intc the roter assen 


device. Maas Wacent=Gura-10n avolds 


”n 
OD 
Qu 
ct 
eye 
(Dp 
4 


Mexirngeeet occurs in the typical clo Tn 
generally avoids the complex mid-tube exchange process 
a 


well as the adverse core-boundary layer interacti 


Oo 
e 
9 
Sa yD 
om 


nn) 


Sebong acceleration force present in the periphery of ¢t 


Boror, Of Clecteic retating mechinery, coupled with th 


w 


heating due to the electrical losses, induces stong loo 


oO 


ereculaction. This heat can then be rejected sither to the 
th 


(0 


shaft of the rotor acting as a heat exchanger itself, <*c 


ge@=eepasSing =hrough the hollew rotor region and the a 


Jo 
ry 


passing through the end-bell region of the rotor ends, cr 


ct 
oO 


_ 0 


borh. The shaft, acting as a heat-exchanger itself, can b 
ecolee with an internal, rotating heat-pipe (or two-phase 
ther mosyrhon) Which could transfer the hsat axially to 
another keat sink outside of the motor/generator casing. As 
ar alternative, the shaft eculd be internally ceocled by 
CmmrconveCcion Of water (cr some other fluid) o an 
Heemeneat Exchanger through th= nen-coupling end of the 
Poeene tse Of an axial rotating seal. The reliability 
and the casualty-control aspects cf the axial rcetating seal 
Make it preferable tc the radial rotating seal mentioned 
earlier. 


Theoretical analysis of the loop-flow has used simple, 


ore-dimensicnal force and energy balances, and need net be 
repeated here. The closed loop is shown in Figure 2.6 as 
conceptually applied to a rotor. Of course, numerous loops 


would be designed into the rotor to carry away the heat 
generated. Inceorporaticn of a one-way valve was suggested 
as being necessary te insure one-directional flow created by 


the mean temperature difference and the pressure differences 
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Figure 2.6 Rotating Loop Thermosyphon for a Rotor. 


Peeeene Lamers of the circuit. Experimental testing was 
Sonmaucted 4n 2 test rig similar to the rig proposed Eor the 
evaluation of the systems studied herein and now under 
Gonstmiction at DINS RDC. A modest range of rotation rates 
was employed in the Merris and Davies test rig, typical of 
most mctcr applicaticns. The details of their test section 
were very similar tc the model used for analysis herein: 
25.4 mm (one inch) diamet2xr copper rod through which a 6.35 
Oe (0.25 inch) hole was bored. The experimental results 


were placed in the form: 
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miceratsc Of angular acceleration to gravitational eacce 


@lera- 
Eeonetac) WaS the authers' way of including the rotational 
etfects (previous authors have used the rotational Grashof 
Number). The final correlation was: 
Nu = 0.1505 Ac0-735 Re2.4%5 PryGr (Ze25)) 
for the range of 103< NuGr/Pr <106. 
Micwrceic»cal ccnsidsration in the analysis of the 
closed-lcop thermosyrhon is the temperature rise in the 


heated limb (conductor bar). The model described in Chapter 
I was analyzed by using the governing relationships from the 


Davies and Morris [Ref. 45] paper iterated against the Wocds 


euamren=:S (Ref. 31] COE rslation. £05 laminar, fully- 
develcrsed flow in a-vheated circular tube rotating about 2 
parallel axis. It can be seen from Figure 2.7 that the heat 
transfer increases with ae cee higher acceleration 
TABLE II 
Temperature Rise in the Closed-Loop Thermosyphon 

REM Temp. Rise( C)* 

300 209.0 

1000 I3%. 5) 

3000 41.0 

3600 55S 

5000 2a e0 - 

* for a cold-side temperature of 45 C. 

errect) . For @he cenditions on the model which constructed 
Figure 2.7, the values of temperature rise is shown in Table 


II for various RPM's. This clearly shows that this method 
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ee cocling, Lf selected, wonld be inappropriate fer the 
Motor device (2.e., lcw RON); ene ese other hard. — 22 could 
Beever y scascnadlo for the generator. DiS See een fac 20> 


was assumed ¢95 be 64/R8e and ne corre 
ie eM 


bh 
42 
a 
ray) 
He 
fu 
mB 
te 
s 
oO 

x W 

\.te 


{ 
(Laninar) are included, the temperature criss 
n i11 be the s 


epe=eqsed slightly. The wtTe W a 
temperature rise at the higher RPM's will still be accep- 
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a 


saple. The actual design of the closed-loop thermosyphon 
Welt rteguirs an optimizetion of the heat-trarnsfer calicula- 
feonseand the fricticn-factor calculations. This optiniza- 
tion may indicats that a larger, single, closed-thermosyphon 
eopmeror a Conductor bar group is preferable over individual 


€eDquc=cr Loops. 


C. ADVANTAGES/DISADVANTAGES OF LIQUID COOLING i 


The main objectives cf cooling are to remove the heat 
created by the electrical iosses and improve the operating 
CeErCLency, extend the lifetime of the insulation, and 
reduce the overall weight and size of the device. The 
feature cf liguid cocling (vs. gas cooling) is «he superior 
capacity of liquids, especially water, to remove heat with a 
small rise in temperatures. The principle disadvantage of 
Ms2ng liquids, or gaseous fluids other than air, is the 
Gesign cf the piping, coupling, and secondary heat exchanger 
required to support the device cooled by the fluid. Since 
ene of the main objectives is to reduce weight and size, the 
fluid chosen must be able to remove enough heat per unit 
mass =O accomplish this goal within the davice. 
Additionally, since the application under evaluation is 


marine cembatants, tke fiuid should be non-toxic, ard <he 
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b 
emeunacceptable choice fer the marine combatant due +c its 
Pm@eeOenve Tevure and since the sizs reduction cbtainsd is 
met sufficient tc warrant system development. Of eke 
momeereng tluids, the obvious choice for its high heat 
emeacity, ave labili+yv, and safety, is water. The use of 
Weter has been tested both in stators and, on a limited 


Scale, in rotors. 


© 


he principle problem to cvercome in the use of water is 


x 
in the rotary sé¢als needed *o direct the fluid inte/ceut o 
the machine. Additionally, when the system inside the 
device is eénalyzed, the ability of a design to withstan 
Gamage is questionable. therd=-siqusewould b= -ceitically 
devendent on the fluid within the rotor for oparation, even 
Mea cacYaicy mode. The supply of ccoling water to a stator 
Bomed §fG= be as critical, since any casualty could be 
repaired easily with shipboard damage-control equipment, and 
emergency cooling of a stator could be accomplished by 
external means due to the thermal conductivity of the stator 


housing. For the rector, this leaves a choice: either not 


& 


tilize the enhanced cooling and suffer the weight and size 


btw 


nalties, or develop a closed system for cocling the rotor. 


yO 
1) 


he closed system for the rotor is not necessarily isolated. 
A ciosed-locp, thermesyphon could be placed in the rotor, 
Beene teseang the heat from the rotor to the rotor shait. 
This locp could have additional heat transfer to air passing 
cmeotwgnecie core of the rotor. Further, the heat trans 
ferred to the shaft could ke removed to the ambient air or 
an external heat exchanger utilizing the non-coupling end of 
the motor Ey the use cf circulating water through the shaft, 
cr by a rotating heat pipe as seen in [Ref. 48]. 
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Cee camceo ince. one t 


GC 9 
peen studiec by numetcus authors, is the use ci 
% 


greatly enhances tte survivability of the device 


damage-ccntrol sense. Bach heat pipe is indspendent of the 
others anc damage to a few, unless all in the sam2 grouno of 
eaqgecent conductor bars, would not jepardize the device; 
indeed, even if all were in the same group, limited opera- 


Tron Gould continue. 


. 





IiI. TWO-PHASE FLUID COOLING 


Brewconcept Of util: zing tne Patent heat of ve 
tc «ransfer heat is widely used in heat exchang 
where, Efrem househcid and industriel boilers and 

heaters to marine prcpulsion units. Air-conditioning equip- 
ment and heat-pump units also employ this idea to a 

MeoomtciyesQtil1ze the lew temperature of vaporizaticn of 
fluore-carben compounds to remove the heat from the ambient 
and release this *thrcugh condensation under pressure, hus 
transferring «he heat in the desired direction at the cost 
of the energy to pump and pressurize the gas to the point of 
condensation. In each of these two-phase applications, the 
transport of the working fluid may require a pump, cr it may 
ke accomplished by the buoyant force that exists in an 
acceleration field, such as that due to gravity. TREES is 
the principle that is employed in the closed two~phase ther- 
mosycvhon. TELS device is simglar to the classical “heat 
Pepe" ane the only difference is that the "pumping”™ method 
Eemege heat pipe uses Capillary action instead of accelera- 
me2oneroress t> pump the fluid. Thus, a "heat pipe" used in 
a rotating reference utilizing the acceleration forces to 
eeanecrort =rs M#guid from the condenser section to the 
evaporasor section (in lieu of a wick) is, technically, a 
two-phase thermosyphcn. Figure 3.1 illustrates the opera-~ 
m2onmer =ne closed two-phase thermosyphon. The evaporator 
end must be in the direction of the acceleration field for 
ene sliguid £o be transported from the condenser. Figure 3.2 
Meliecrates the operation of the heat pipe, which utilizes 
PuewcGapot lary action of the wicking material to transpert 
the liguid from the condenser to the evaporator. The 
remainder of the chapter shall use these definitions fora 


discussion cf both devices. 


We 





Mecelenation Field 





Figure 3.1 Typical Two-Phase Thermosyphon. 


Acceleration Field Not Required 


\ 





oe 
\\ Liguid 


Flow 


Figure 3.2 Typical Heat Pipe. 


A. HISTORICAL DEVELCEHRENT 


The development cf the heat pipe (which will apply 
equally to the develcpment required for the two-phase ther- 
mosyrhon) was begun by Gaugler [{Ref. 49] in 1944 and by 
imeceetnen {[{Ref. 50] in 1962. Trefethen was working on 
cooling methods for spacecraft (zero-g environment) 204 


General Electric Cempany and discovered that capillary 


aS 





pumDding can be a very valuable area for further develcoment. 
Working independently, Grover, et al. (Rete ol jee eo Ver cee 
Pee cor ecen-; he widely published the first applicaticns c* 
mee d¢vice ard gave it its name. imeus, et ai. Hee iiecg  oae 7h 
have rerorted thi historical development and have given 
Sample calculations for the design of a heat pipe. Cai 
(Ref. 53] also reviews the theoretical developnen+ of heat 
Bipes, and frowides the design procedures, including seaple 


Paetewidt2 cns, TOr their utilization. 


Be. CONVENTIONAL, CLOSED TWO-PHASE THERMOSYPHON ANALYSIS 


AS previously mentioned, the appealing facters associ- 
Beedwew2th =he two-phase utilization include: high heat 
fluxes asscciated with phase changes, lower ‘temperature 
gradien+s associated with these processes, and reduced 
weight of the two-phase system over the single-phase liquid 
system. 

' Research by Cohen and Bayley [Ref. 54] (referred to and 
Giscussed by Japikse {Ref. 44] ) found that the amount of 
liquid filling the system functionally affected to the heat 
Pawmenen 2m both rotating and static tests. They found that 
heat transfer increased as the percent liquid in the systan 
was increased +o approximately 1.5% by volume, then 
decreased +o an intermediate value, and finally increased to 
approximately the same value as in the 1.5% case. This was 
related to the following process: eee ne Vow Filling situ- 
aclon, cordensate returns +0 the evaporator section and 
Por@e doenin Ei lmon the walls. It is in this film that the 
heat transfer occurs and the phase change takes place. ir 
the case of the completely-filled evaporator, a liquid pcol 
region develops, and a condition of nucleate boiling exists. 
Pemececrlting 1S insufficient, with regard to the heat flux 


bBeang trensrorted, the pool/film will not continue to the 
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Pvemcte=uyc heated section and "dry-out* will occur. is 
Wary-cut" will assume a drop-wise instead of film-wise chass 
change. Hue geenetry Of The regeons in which ¢: 
dzkely tc form, influences the results. 

Les and Mital (Ref. 55] conducted experiments with an 
electrically-heated, water-cooled thermosyphon using water 
myeomerrcOnu) eS th] test tluids and varied the filling quan- 


and heat r1 


G 


XoCG yk SSserhe 
Sane ( ie a Ss 


ren (o2T 
£ the heated section to the cooled 


ssu 
he length o 

Om). Mae besulc of the filling quantity or heat 
nefer was the same as that reported by Cohen and Bayley: 

@sing heat transfer with filling to a point and then 
decreasing beyond that value to an intermediate value and 
increasing to the case of the completely-filled evaporator. 


Bewece-cc- Of decreasing L./l, was to increase the heat 


een c=sr Withir the range 0.8 to 2.0 and the advantage of 
larger ccndenser area was evident. The heat-transter co¢ti- 


E2cten= Was found to increase sigqnificantly with increasing 


S 
mean pressure due to at least three factors: 


1) since the mass flcw for a given heat flux is nearly 
Gonestene and density of the liquid increases with 
Dressure, a lower pressure drop (and lower AT) is 


necessary fer the same heat flux; 


Ze £02 a larger P varies much more rapidly with 


D 
sat’ “sat 


© : for the fluids considered, hence requiring smaller 
sa 


AT's at higher pressures; and 
5) for lower pressure drops, a more favorable force balance 
exists on the condensate film permitting a faster liquid 


return. 


Water was ‘found to give heat fluxes superior to those of 
Seon, FOr for the sane AT, due to the larger values of 


latent heat of vaporation, and thermal conductivity. 
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where Glew t7G8-) Wyo ye +9378 y* — 172 y* In y , 

even y (17/2 in ye -1/2) + 1/8 y* #+y2(t/2 ln y -i74) =178 
and y = 1- O/SR. The simplifying assumptions (forces due 
tc vaper shear and momentum changes were neglected in ths 
force balance) causé an error (of as much as a factor of 2) 
Zooccur, but qualitative behavior is correct as far as 
Li/lLoe atc working fluid are concerned. This develop- 
ment is also included in the survey by Japikse [Ref. 44]. 
This heat transfer is graphically presented in the Lee and 
Mital paper in their Figure 12, shown here as Figure 3.3. 
This clearly shows that the quantity of heat that can be 
removed by a two phase water device is large, even at a low 
saturaticn temperature. [PNescepescral sor ind*cates that for 
Milestone —Inereasing ratio of heated length to cooled 
length, L+, decreases the maximum heat trensfer. The heat 


transfer is also very sensitive to the operating pressur2; 
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Paeguce 3-3 CoOmpariccn © the Experimental Results+with the 
Analytical Prediction-from Lee & Mital paper. 


the heat-trarnsfer coefficient increases appreciably with the 
méan cperating pressure. The temperature drop increases 


with the pressure drep aleng the length of the thermosyphon 
ration pressure gradi¢ent). As the ovres- 


tube (and the satu 
sure increases, the vapor specific volume decreases, 
resulting in a decrease in the pressure drcp; the mass flow 
Sete of fhe vapor is essentially constant for transferring a 
gGawen heat input vette. In ‘the previous papers, the 
ver*~ically-oriented +thermosyphon was considered and the 
EVErent model (for scme of the orientatiors) requires the 
device teyb= oriented in he horizontal direction (perpendi- 
Gulez to the applied acceleration field). Aithough the Lee 
and Mital paper was considering a vertically-oriented, +two- 
phase thermosyphon, the general results ares th2 same for the 
HOr2zencal tube. 


el 








oe LISITS OF OPERATICN 
Wien CORSiGerzing the d=sign of a tvo-phase thermosyphenr, 

Gamo wei-ceasbipe, lindtations to the heat transfer must bd: 
considered; four of thase limits are common to both devices 
(source, encrainment, toiling, and condensing). A qualita- 
tive ccmpatison of these heat-transfer limitations as a 
eypciica Of “he ssture+ion temoerature is givan in [ Ref. 48] 
aoe avs Shown in Figure 3.4, 

< BOILING 

- LIMIT 

< 

c 

Ge 

had 

\ 

4 ENTRAINMENT 

< LIMIT 

ze 

-_ 
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~ CONDENSING 

xa LIMIT 

SATURATION TEMPERATURE 
Pigure 3.4 Operating Limits of Rotating Heat Pipes. 

The sonic timit andthe boiling limit have been readily 


analyzed. 
not 
these phencmena 


research is needed. 


The entrainment and 
as well-known and little, 


exist. 


the condensation limits are 


tivany, literature describing 


Dais 2S Wan area in which further 


the vavor flow in a two-phase closed 


HR esonve svelocity at the operating 


To 





peéssur2> (also known as the "choking" limit). The sonic 


tame 2S Sepresentec ry: 
Q 285. =5, 44h. . (6523 
max rg ¥ VS ig 
where vis they sonic velocity of the vapor. Lies ene ind 
velocity may be determined experimentally, computed by the 
Betaeionship V_={((3£730) Malice or approximated bv the 
iS * 
Pertect-gas relationship Cetra t/ 2. 
2) Entrsinmert jlimit; the interfacial shear between the 


H2@uid and he vapor will hcld the liquid (which is 

direction) back and starve the 
uid S$ counter-current flow, when the 
relative velocity is large, causes the interface toe beccm2 
unstable which results in waves at the interface. As the 
Mepo= Velecity continues to grow, dzoplets of “J2quid are 
formed at the liquid surface as the shear force exceeds the 
suriace-tension force. The formation of these droplets and 
their suksecuent entrainment in the vapor stream causes the 
PamceciwoL tOta! Stoppage of the flow (dry-out). Ease 
phenomenon is generally governed by the Weber number (the 
Bertie CE the inertial force to the liquid surface-tension 
BOEC Ee )-< The Froude¢ number is also used *o characterize the 
phenomenon cf the drep-wise entrainment of the liquid in the 
Wapor strean. In the application considered herein, the 
formation cf the waves is considared unlikely due <‘*o the 
extreme acceleration field present. Since the flow is 
coun*¢r-curztent and is in the presence of a high accelera- 
tion fisld the entrainment limitation thus will reduce to 
fhe "hold=up"  dimit. Tetwilil remain stratified until 
hold-ut OCCUES, resultin in evaporator dry-oat. 
Experimental determination of the exact correlation for the 


hOld= JO, OL entca naent, limit is requized for a horizontal 


qg 





Sug-m2Oece1nd around a parallel axis, thereby creating é4 


perpendicular acceleration field many times greater than 
Geese. Tae entreinmert limit is gensrally given by: 
zh ty 
Q max Reg yA re (3.4) 


where V_ is when the Weber number equals unity, the Froud? 


P 
— 
S 


mUmDer Scuals unity, or is €xperimentally found. In ancther 
Beppe Gece, Jester ard Kosxy [Reft. 56] used the ratio of the 
eemalechear force of the vapor flow on the liquid surface to 
the gravitational (accelertional) body force upon the liquid 
ige2i:ned as ‘F) Sr olae= tO arrive at an appropriate 
"hold-up" criterion. They correlated 2xperimental data as a 
function of this "F" value, as the flow they were studying 
tEansiticnei from stratified to annular. The result of 
Treis €xcte=>smenm=S was that the flow wes stratified for "F" 
values below five (5) and was fully-annular above twenty- 
mune (29). The value of "F" set agquali to five (5) was 
Perot hkec as Lic €fLitea=.On Lor the transition from streti- 
fied tc annular for the co-current flow case. The seme type 
of analogy was develcred by Collier and Wallis [Ref. 57] who 
Herearcecetne inertia force end the acceleration force to 
scale stratification effects where their criterion was given 
es) WF (volumetric Elux) equal to 0.25. Both of these refer- 


ences point out the need for further rescarch in this area. 


3) Poiling limit: Again, as discussed in the first 
Ee@euGnwe: this Chanter “hn the general discussion on two- 
phase thermosyphons, the boiling limitation in the current 
application is due tc the creation of profuse nucleation at 
the ev2pcrator section and a resultant vapor film on the 
evaporator surface which insulates ths evaporator well, 
Pov noun tay-Olt and overheating. Tie boiling Limit is 
Given by the Zuber-Kutateladze prediction as discussed in 


(Bef. 48] and is given as: 
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4) 

dependent on many haat 
Weztking fluid, the o 

the acceleration fie 
crientation of the vessel, for the model considered, is that 
cf ahorizental tute (perpendicular to the accelaration 
mesa) With internal conden sation. The condensing limis is 
the capacity to condense the vapor on the inside surface of 
the condenser section, and this is shown in Figure 3.4. No 
analysis on the outer surface is presented in *this thesis. 
mcwedd,eth= condenser init is based on the capacity of the 
inside wall exposed to condensation, using an arbitrary, but 
reasonable, wall temperature. (Generally, «the largest 
themrmal resistance cccurs between the ambient and the 
condenser cuts surface. Therefor2, heat transfer is 
limited by this resistance unless the outside surface area 
is substantially increased by the use of fins.) . The 
condensing limit is aiven by an equation from [Ref. 58] by 
Gelecry2Or a hozizontal tube with internal condensation in 
a gravitational field as: 


Q =A ALE { (6 


= h k 3 
max Vv Cc lee aes p19 fqaue Le 


Puse( oe s, yyy se (3.6) 
Foro the rotating referenc2 being considered, «he "g" term in 
equation (3.6) is replaced by the centripetal acceleration 


G2 eR. ivomececctere s “ Inmenis equation, allows for <ths 


fact that the rate of condensation on the statified layer of 
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MeqiscelLe Reg. icable. Phe@eyalue of "“F " dénend= on zh- 
angie »¢ ana&d is tabulated in {Ref. 58}. Its value ranges 
fom Zerc, When She tubs 2s f£u22 (no su 
Seidenmratiton to OCCUL), Oo 0.725 , when the <ubs empty. The 


eagle © 2S shown in Figure 3.5. 





Figure 3.5 Laminar Condensation within a Horizontal 
. Tube (from Collier). 


1] 


sation of Heat-Transfer Limit 


Me calcula<icn of these limics is an important step 
in the design process. 

Bheecalculaticus Eegquzred to anbalyz= the configura- 
t20n cf the current medal are not all weil known. rac 
following appear to be the most applicable to the model 


being evaluated: 


e=—=—-The sonic Limit calculation is from basic principles as 
shower by equat:on (3.3) by ellowing the vapor velocity to 


increase to the Mach value of unity. 


|“---The boiling limit is given by the Zuber-Kutateladze 
predicticn as discussed in [Ref. 48] and is given by equa- 
meeon (3.5) 
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Meme cone ulS reveals no correlation yet for the entrain- 
t in the heated, Nhepezoncal cube rotating ebcut 4a 
GxeouwvenwCCUNnTEr-GCUrrent tlow. Since the ma jority 
es 


previously conducted for this type of flow have 


SOW cue. <8 CDtreéegmment linit is typically below “th? 
Eemtc sand aebovw> the condensirg limit, the correlations 
previously mentioned have been plotted against the amount of 
meme that 2s Charged in =he tube (percentage fill) for the 
three reasonably known values (sonic, DoLit, and 
condensing linits) Loewe Sa GUTS Se 7s Also shown is the 
entrainment limit ccemputed from the Jaster ard Kosky 


feet. 56) Correlation developed for co-current flow. 

homves Dreviously noted, the condensing limit is the 
Bteecr ttat as the mest d2miting wn the design of the ther- 
mcesyphon. The values of OTe =e increased to a maximun 
Bassa on the length of the condenser, SNe) percentage E£ilt 
and the AT availatle between *he constant temperature 
(assumed) wall and the saturation temperature of the vapor 
condenmeing on it. As already mentioned, the resistance to 
heeat-transtfer afforded by the condenser outer walls to the 
ambient may b? quite significant. The ability to design an 
adequate fin arrangement and supply cooling air in suffi- 
cient quantities, whil2 minimizing the windage losses, 
itself is a substantial problen. From the data obtained in 
his rough model, there are some self-evident points to be 


considered: 


----The length of the condenser, ‘from a heat-transfer point 


cf view, determines the heat-removal potential, 


=~ wewamemre Of 1iguid fill can be optimized, 
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For the model being evaluated, and neglecting (as 
the Jimitino value) ehewenepatnmer: limit, the closed, 


tea-phas¢ thermosyohen, operated at a saturation 


a 


tamperature 


wall temperature of S0°C with rotation at a 0.381 m radius 


at, or above, 100°C and with the condenser . 


tia 
(Dp 


Qu 
oo 
cr 
fo 


chien 


‘3 
co 


and 3,600 RPM can transfer the required 50 W of power with 


the condensing limit being the limiting value. Prqus< 3.7 
15 based on a heat load of 50 W to be equally divided 


between condenser sections in both end-bell regions (i.6¢., 

approximately 25 W each per condenser). Fos the 3600 RP 
model the percentage fill may increase to approximately 40%, 
and still maintain the required heat flux. The effect that 
the REM has with regard to the temperature rise in the 
heated section, makes the chcice of the clesed-loop *herno- 
ByOnNOn, 2COr Or Ehe motor applicetion, but quite suitable 
for the generator. Further Dackground on the theory of 
two-phase flow is available in publications by Wallis 
fRet. 59] and by Collier & Wallis in [Ref. 57}. The f¢easi- 
bility of the design has been shown by Corman and McLaughlin 
fmet. GO) and by Groll, et al. (Ref. 61} through experi- 
ments with heat-pipe cooling of electric motors in the 
scheme primarily reccmmended, although the difficulty of the 
liquid transport at lew retation rates was noted. A diagran 
of the cooling scheme such as recommended and tested by 
Corman & McLaughlin and Groll, et al., is shown in Figure 
S26 (Eeproduced trom [Ref. 60} ). 


84 





"suayos HbutToo) sdtq—-zeseH TeotddAy 9°€ SInbty 





SID 7 VEIN 
nee = ce ~ 






10104 











= mat 


= i 
WAL. EZ NN 
o i 8 


























85 





. ee ene: aseyd-OaL 
PSssoTD 8Y4z UT 


ns 


~J! ,U! 





aywyi Buyeuspuo7 


cE “e = sd 


q) 2OrF*O = ENTPOY 


COSE = Wd 


aywy Bur prog 


Sica ee 2a Apwry] DyUOG 


—. aes 
- 





q1ut7q Buyuyooyug 


ZTUTT AdJsuezry_—zesy Lee sanbty 


ey 


or 


at 


ont 


TBQoL 


CD) 422 


86 





As discussed in the paper by Green2, 2t al. [ReEE. 6245 
Esc nest foe el Vv 

and weight is well k 

advantages of this 
constraint “io be overc 
cf these systems is ther 
devices. Ine t976, We Ds 
increasing the load er unit weight 
demanded that the inherent heat created byt 
losses be removed via improved cooling technigue 
also elaborated on «he "complex arrangement" of th 
and asscciated equipment required to accomplis = 
Bec eased cooling. Further, he stated that "There i 
currently STsumeietene technical infotmaction aveilabk ie 
Scrmeco Nak eG  CONtaden:  predictcons of the effect of 
Coriolis and centripetal acceleration on flow resistance and 
fewer SanStes 9 17 the complex “coolant cizrcuits...". Sines 
then, Merris and his co-workers have attempted to fill the 
literature gap in this area, and have uncovered a great many 
acts tegazding this flow structure; they have also uncov- 
ered a great many guestions for further research. 

AS shewn as early as 1955 in @n article by Sir Claude 
Caubepictemouy,. the principle factor in the failure cf two 
very expensive generating sets was the inability of the 
devices to dissipate their cwn heat. Im the late 60%, 
Feach [Ref. 65] wrote of the successes that th English 
flectric Company had with regard to liquid cooling of the 
large generator rotors being used in their 500 MW devices 
and predicted that the devices would reach 1,000 MW with 


mis technology. This was followed in 1970 with a second 


og 





Sie OS} eNne=e he aujlined the ercv=h of 


etmtien whe ligquig-cocling 


€ 


R 
izes 
. 


and c 
Sait Sedecjmone NSat fLitars Ye susocr= 
is 


she devices ing planned in thi country which excezeded 
TOO) Ms The losses developed in the large-scale, high- 
power, superconductive devices, currently heing developed 


for large-scale implementation, require extensive use of 
cooling at cypogenic temperatures in o 
the sucver-ccenductive stats in which they 
Sion of this applicaticn is presented by Schwartz and Foner 
(Ref. 67], in which a number of the schemes discussed herein 
are shown. 
The previous chapters have discussed the ability of 
ertnral, Forced-convective cooling =O rencve the 
internally-developed heat loads in ealectrical, rotating 
machinery. The reduction in size, the increased efficiency, 


and the prolonged life cycle are all very important reasons 


to utilize advanced-ccoling cencepts. The use of turbulen+ 
fiow is feasible at all FPPM's; SnererOnrs, both motor and 


generator applications would benefit from forced-convective 
Ba iea Cocling utilizing turbulent flow. Laminaz-flow heat 
transfer is feasible only at higher RPM's, which tends to 
ia@eatethe appiicaticn of this method in regard to motors. 
Piewacetel transition from laminar to *urbulent flow is not 
well known and ‘further investigation is required to dster- 
Sense che tzansition foint; eNescEensit=ton may occur at 4 
flow and REM acceptakle for the motor application. However, 
cere utzlizaszon cf external, forced-convective cooling has 


some major difficulties associated with it: 


1) The source of the external cooling fluid must be 
easily accessible, dependable (in a casualty-control sense), 
safe, volumetrically small, light, and should be reasonably 


inexpensive. 
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eee fa atseli should be of high specific heat, 
tz2q37e, inexpensive; and more importantly, NeON=TOXKLC and 


By Fhe gwetkhoad of transp 


c 
through an #xternal heat exch 
= 





back, Should be extremely Siant to 
external damage. perce pemmman 4d OL Sev lous conside2a7icn 
weeeemo ne sels ebliiaty of conventional, rotating steals is 
musstecnable in the least. 

wre alternative to the use cf external, forced- 


convective cooling is the utilization of a closed device: 
(a) the clesed, twe-phase +thermosyphon, (bd) the closed, 
rotating-locp thermosyphon, (c) the closed rotatina-loop, 
two-phase «hermosyphcn, or (d) the heat pipe. These alter- 
Natives remceve the necessity of using radial rotating seals 
in the design of the cooling systen. 

The closed, two-chase thermosyphon, as shown in Figure 
Eeopmemnlsct Neve 2 sufficient forced-air circulation through 
*he end-rkell regions to remove the heat... This will require 
cptimizing the length of the condenser sections, the gsom- 
eemy or the fins, “the Bamber of fins and the direction of 
air flow within the device to maximize the heat transfer 
while minimizing the windage losses. The closed two-phase? 
thermcesyphon is feasible at the higher RPM's that the gener- 
ator operates at; additionally, this method could also be 
feasibtle in the motcr application if a heat pipe (i.2., 2 
thermosyrhon whose walls are lined with @ capillary wicking 
material) is used in lieu of the two-phase thermosyphon for 
che lower REM s. 

The closed, rotating-loop thermosyphon requires an addi- 
tional heat *éxchange through the shaft of the dsvice, 
possikly supplemented with forced-air convection through the 


hollew retor region and by the exposed radial secticns of 
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, Closed ‘Water spr s = 
may also te cooled by for ec a 

g, with axial, packing-type seals, axter 
@evBce as skown in Figure 4.3. hes. 1. p 


3 * am } 3 1 a= > a - F aes =, 
pe@rme Could USe elther Single-phase or f£Wwo-phats cperation. 


The closed two-vhase thermosyphon could also be extended 
along che shaft ‘threugh a bearing plate arrangement and 
Poemangeet=. Could Ee supplied external to the casing, 
S#egplifyina the internal configuration, as shown in Figure 
G4. 


The ability of these devices to cool effectively has 
been proven and their usage in this applicaticn is theoreti- 
cally feasible. TiesexeacesCOrbelations for the calculation 
of the performance values for these devices, especially ina 
BOkatang rererence, will ftegquire further research. 

The analysis herein hes shewn the ability of th> cocling 
schemes discussed to remove the requisite heat load in erder 
eemexrcend the Jite cf the instillation, increase the ~ffi- 
ciency, decreas¢ the weight and size of motors and gen2ra- 
fOrs . This study graphically presents the applicable 


Cepcela=son]= for the determination of the heat transfer 


if 


mectinea ctOteting, electric device. It alsc discusses the 
linitations of these correlations. It has been shown that 
closed, two-phase systems of an inherently higher reli- 


ability are feasible, and deserve further evaluation. 
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¥. CONCIUSEONS AND RECOMMENDATIONS 


ee ee eee eee ee ane ee te So we ee ee we ee eee re ere a So ee eee 


BR COCNCIOGSIONS 


M 
A» 
G 


ras pon the information reviewed and analyzed in this 


Baosts, the £0llowinc conclusions are mads: 


WO =cea-—convective liquid cooling is a 
BeeeesOr I0proverG “~he EEiticiency, reducing the siz d 
weight, and extending the life of high-power motors and 
generators. The sein disadvantage of forced-convectiv: 


b2qu2¢d ceolzng is the requirement of radial rotating seais. 


2) Closed-loop liquid thermosyphons are feasible at 
hign REN. They offer improved reliability over forced- 
Gomvect.ve liquid ccoling which utilizes rotating seals. 
However, this method requires the us2 the shaft és a 

econdary heat exchanger 


Bj the closed, two-phase thermosyphon is feasible at 
mee higher RPMs of the 2nerators. This method requires 
forced-air convection through the end-bell regions or the 
extension cf the devices aleng th2 shaft for external 
forced-aizr convection, both of which will add to the windage 


losses 


ive cre low RPM applications, such as motors, the use of 
a heat vipe Wiel epe Tequired to overcome the inability of 
Peewleumaccclcoration field to transport the cooling fluid. 
Both the closed, twe-phase thermosyphon and the heat-pipe 
require forced convection at the condenser ends +o remove 
the heat either internal, oor external to the device. This 


may add to «he windace losses. 
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sole a2 the litertature fer: 
elec eeeGOUNcer-CUl= ens Grtrainment limit for tvo- 
phase thermcsyphcens under rotation a 
Dyetn= transition from 1 = LCP cus ou Lem 2 OW 


iene = Otating*systems. 


E. RECCHAENDATIONS 


meer ollOving Crecemmendations to ~continu> this work are 


macd2: 


imtescinG “Ore the f£oreced-convective liquid cooling 
scheme should be acccmplished using “the DTNSRDC test rig to 
confirm the analysis and the validity of the available 
@rcelations, GE tom ODtain PEOpSE, correlations for the 


specific gecmetry and utilization being projected. 


2) An experimental program should be devised and 
conducted tc determine the entrainment limitation in a hori- 
zontal, heated, channel rotated about a parallel axis with 
Come ece-curcent flow cf liquid and vapor. 


3) An experimental program should be developed to 
determine ‘the optimum condenser length and proper fin 
arrangment for the use of both heat pipe and closed two- 
phase thermosyphon systems for rotor cooling. Consideration 
Compe nceructural integrity must be included for the 


condenser end-lengths. 


YW) An experimental program should be developed for 
correlating the laminar-to-turbulent transition during 
Beueed=couveet: Ve 1iqguid cooling in a heated, horizontal, 


Eoaae 2d p2De. 


5) An experimental oprogram should be developed to 


- 


modify tke inside gecmetry of the closed-loop +hermosyphon 
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wt 


a 


es 
14 
tN 

vd 
cr 


(emexae ll, Closed, <WO-vhase “hermosyphor} <*0 op 
a 


eyic 


fs 


Neo tee =rictior losses wathin the 


w 
o 


6) A dstailed analysis should be made of the shaft 
cooling fotential fer the rotating-loop thermosvphon with 
regard to the heat-transf enhancement cf the shaft to the 


nsifec 
fluid and the mechanical desian of the sheft. 


7) An e@xperimental program should be developed to 
a@palyze the clossd, two-phase thermosyphon with radial 
a 


1 
Llons and the closed-loop two-phase thermosy- 
+wo devices may 


experience start-up balance 
vroblems. 


8) Various configurations should be analyzed for their 
petential for heat transfer to include a combination of the 


h 
systems described herein, as depicted in Figures 4.1 through 
U.4u,. 


Syme tsewcencep: cf Noditying the gqaometry of the rotor 
Bandangs wechin these devices in otder to mininize heat 
eomamiceionepath lengths to the casing, shaft, stc., “and 
thereby the internal temperature gradient, should also be 


investigated. 
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CONVECTION IN UNIFORMLY 
BREATED HORIZONTAL PIPES AT LOW RAYLEIGH NUMBERS 


a 


The following develop 


_ 


ment, BYeenmor On [ Ret. 17 }s is 
presented kerein fcr completeness. The steady laminar 
Mer2cn Cit a fluid in @ horizcntal circular pine of radius a, 
the walls of which are heated uniformly so that a constant 
temperature gradient +t is maintained in the directicn of 
ceeomaxec, ths flow. wili be referred =o in cylindrical cocr- 
dinates (r,6 ,Z) with ¢ measured from the upward vertical 
and the z-axis along the axis of the pipe; the velocity 
components are denoted by u,v, and w. The effects of dissi- 
pation and cf the pressures term in the energy equation will 
ke neglected and that variations in the density due to the 
temperature differences are so small that they only effect 
ene buovarcvy t¢rn. Thermal conductivity and kinematic 
viscosity are assumed to be constant (Which will 2neroduce 
G@igmeitae-ive Grrors into the solution, but shouid not change 
Peeeqenetal character). 


The equation of ccntinuity is then: 
dru/fdor + av/ogd + Orw/oz = 0, (A.1) 
end cre energy equaticn is: 
Meus) VY/El(oT/oo) + w(dT/oz) = k Vet, (A.2) 
where V2=o2/gr2 +1/r(d/ar) +1/62(d92/o¢?) +392/az2 and T is 


the local fluid temperature. The momentum equations can be 


wrest en in 2ne form cf equations (A.3) through {A.5). 
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U(du/or) +4¥/r ( dU/ 50) tw ( du/9zZ) “v2 /c = -1/p(dp/dz) 
ey ay =< —27 r= (5%/ 30) — G(T -T) cosy : (A773) 


u(3v/dr) +V/ro (dV/9) tw(ovV/ dz) zuv/o = -1/pr (dp/ do) 


pate 2275S 40 / sa -¥ /> 2+ See ola) ; (A.4) 
U(dw/dr) +V/0 (3W/390) tw(dw/9z) = 
-1/0(9p/9z)+v yew. (A.5) 
In equations (A. 3) and (A.4), the buoyéency force has been 


calculated relative to the fluid at the same level adjacent 
*o the pipe wall and the remaining distribution of force has 
been absorbed into the pressurs3, p. 

For steady convection sufficiently far from the pipe 
@xemirag £0 avoid inlet-length effects, the temperature 
throughou* the pipe increases uniformly with the distance 
along the pip2 axis. Hence, the distribution of buoyancy 
Poeee Ah sections of the pipe is independent of Zz; as the 
secondary flow is caused by the buoyancy forces, the flow 
field must also be independent of z. It follows that there 
should be a similarity solution with u,v,w, and T.7T as 
functions of r and » only, and with p=yYz+P(r,? ) (where P 
contains the «erms aksorbed into p). 


The continuity equation reduces to the form: 
g(tu)/ar + ov/ao = 0; 


hence, a @4imensionless Stokes stream function yw can be 


antroduced in such a way that: 


ru/v =ob / doe V/ Vee = 07 Obs 
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PicetetBl=garion=s (h.1) threugh (A.5) aay be reduced to 
qd 


Pesos e-nercns! form by the transformations r=aR, w={v/a7k, 
ead Toc TstaPr e Where Pr is the Prandtl rumber {v/n}. age 
the pressure is eliminated bstween equations {A.3) and 
(A.4), the resulting momentum equations are: 
my U +178 (0079 (9/04) ~ot/ 96 (3792) } vey = 
Ban(sao7orjsing + 1/8 (d6/dd) COS 0}, (A .6) 
VER I1/R{ (ov/dR) (3/09) - 
(9U/90) (O/AR)}W +4Re =0, (A.7) 
ani where, 
ieee G=7an= tI Rofo) te R= (eae 7 one) e 


The Rayleigh number, Ra, is given by ggrat*/(av ) and the 
Reynelds number is the normal Reynolds number, Re, for 
Foiseuille pipe flow ktased on the pips diametsar and the mean 
yelocity across a pipe section is (a/v) ( ya2/4y ). 

The solution of the energy equation depends on the tenp- 
erature cn the boundary condition at the wall of the pipe 
Pipes will usually have reasonably-thick walls of material 
aren heb Wal Conductivity much higher than that of the 
Pad", so there will be little variation in wall tempera- 
ture around the pipe circumference. This will be specially 
so for slow rates of heating when asymmetries of the flow 
will have small amplitude. Hence, writs T =T0+ z, where TO 
Powethomewoll) temperature in the section containing the 
Caigan. Although there is uniform heat transfer per unit 


| ol 


ength of pipe, the local heat transfer will be slightly 


Q 


reater near the bottcm of the pipa than near the top. With 


c+ 


his assumption, equation (A. 2) z~educes co the non- 


Ou 


imensional form of equation (A.8). 
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V2g + (S/R) ( (907 OR) (99794) - (90709) (9879R) } +H=0 (A.8) 


Weve e, ana 6 Gre Zero on r=a; 
Up eeonGgs ale tanite on r=0. 


{A.9) 


° 


O62 Of Scaquaticons (A~6) zheeugh (A. 8), Setis= 

PyeesGg CopGil-~10nS Of ({(A.9) is considerably difficult, succes- 
n e solution can be obtained by 
expanding wb, W, and 8 as power series in the Rayleigh 
rumber, provided that this is numerically small. Supposing 
that: 
WRay +tRazy tose, 
W=W + RaW ,+Ra 2Wit..., 
9=0,+Rad,+Ra20o+... « 
(A. 10) 


The leading term ¥, c£ wp must vanish because there is no 
@erculaticn when A is zero, but an is not zero since the 
difference in temperature rFetween a fluid element and the 
Merqutkotang wall is rropertional to t6. When relations 
Wenewe) are substituted in equations (A.6), (A.7), and (A.8), 
three sets of aquations for the funtions v5 sw if and oF 
are oktained by equating coefficients of powers of Ra. 


From equation (A.7), the basic equation is: 
v2, + 4Re = 0, (A. 11) 
Woich has a solution: 


#-= Re(1-r2) (A. 12) 


Bamenecatice22s the conditions W,=0 at R=1, W, finite at 
R=0. This is ordinary Poiseulle flow under a pressure 
gradient -4oFRe/ar in an unheated pine. 
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Becomes sce Of he equeticns derived from equation {A.8}) 


for the temperature distribution is: 


cD 


eeemec- =O 25=0 at R=1, ©, finite at R=0, and the solution 
is axially symmetrical as no account has been taken of 
G@eavicy a= =cnis stage of the approximatior. Mifs 2S Sdtis- 
Bred by: 


pa C7 16) Bet l—R-) ( 2=Re), (A. 14) 
Meee. as the customary solution for forced convection. 
The first-order approximation for w satisfies the equa- 
enOT) 

V4 =(399/9R) sin ¢ (A. 15) 
obtained as the coefficient of Ra from equation A.6, and the 
boundary conditions 93y,/9R8,9V,/9¢=0 on R=1, and R-!(9V, 729), 
mayo R relain finite at R=0. If ¥; is assumed to have the 
tommy, (kK)sin® , the dependence of equation A.15 on ® is 
eliminated, and WR) is easily found: 


W, =~ (1/4608) ReR (1-2) 2(10-R2) sing . (A. 16) 


umereorly, | tn fiLrst-ordsr approximation to W,satisfies the 


equation: 
Ge 7) (90,790) (04 ,/9R), (A. 17) 
Seance >sOurdary COnditions, W,=0 at R=1 and W, finite? a* 


R=0. The dependence cf equation A.17 on > is eliminated by 
raka 1G W.=W, @)cos¢ , whence 
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W=- (1/184320) Re2R (1-R2) (49-5122419R4-R6) cosh. (2.18) 


Miemerfe--OLSer aprproxiaMecien to 6 satisfies *+h= Ecua- 
Z2O0 : 
v269 = (0/8) (9), 790) (39 /dR) -W (A. 19) 
Bere TOU Ca = y conci*tions, Cees eee = beams Lu ee oe 
1 
R=0. Hence, 


9, =- (1722118400) Re 2R(1-R2) ( (381413255) -(354 416755) R2+ 


emiguote2 5G 7) R*=(2942000) R&+( 1-25) R®}ccsm . (A. 20) 


This  ces@pletes the approximate solution to the first 
Pecos in Ra fOr convection in a horizontal pipe, which is 
heated so that there is a constant temperature gradient 


weong =<ne Wells and uniform temperature around the girth. 


mt may te roted frem equations (A.10), (Ac t2) 4 (A.14), 
feeloeta.c), and (A.20) that the full convection solution 
depends ¢ssentially ¢n the product RaRe. This is clear if 


ineeet= crecatled thet Ra is proporticnal to the increase in 
wall temperature along the pipe and R2 is proportional toa 
G@@atecteric<i.c velocity cf the flow along “*“he pipe. An 
increase in Ra means that the fluid will be carried through 
a laraer temperature differences, and, hence, there will be 
increased buovancy ferces; but the same effect can be 
produced by an increase in Re, and hence in the velocity of 
-he Wain flew. 

The two most interesting consequences of such a solution 
are the modification of the heat transfer due to secondary 
flow, and the effect of haating on the tangential stress at 
Ene Wal! (c= on the rate of flow for a given pressure 


gradient). However, as both W, and oF vary with cos, 
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Meeens= the flux ner the heat transfer across the whole 
Eeceten Of the oipe will be Bengsa PlDve whe  —2tss=-crrst 


mee San 


wee = (178) (5 ( VE ot )ZO(R oO )} 
Tce yo R) Singt (17R) (36739) COS ’ 
24 =(1/B) BCH V9 (Re OIF CIZR) (804 VI7IRL OES 


T25_=(G/B) (3(9, oh) /3 (Reo) }+ (O/R) (318 oY, /3 (RO) HW, 
eZ: a 
(AS 1)) 


The enly difficulty in solving these equations is numerical 
tediom and che remaining second-order solutions are 
presented in (Ref. 11}. 
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SIMGLE-PAASE ANALYSIS 


The following prcearam was w 
Tne vVerious correlations 


forced scsiraclsa-satasse onvection in tha rotatin rota 
ae reread Wy Sn ese eet) oS Vere Se i) a re Che a Cle Een 


disscussed 


Bodification to 


az was 


Sdb-Erograms for the 


written 


In Chapcer 


OMPOUTER PROGRAM AND RESULTS) 


Srewen £O give the rEresul<s of 
EOme soe heat transtsr for the 


Tok. The program is so written that 
Peelide Gene: “correlations iS) quite easy. 
See —SestGun fOr erie "HP=9326 “Conputer. 
graphics are not shown. Sample runs 


are included immediately following the listing. 


-Pages 


-~Fages 


=Page 


105-109 
110-111 


2 


Program Listing 
Thermophysical Properties (Commcen to 
Appendices B through D). 


Sample Calculation 
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‘10 


aes 
oCecnocoe 


co) 


i) 
P) 


= ‘ = 
- AME 


ne, 


el ww a vel 
Cerrone ott cevae enna 
Brome NEO api me Ariat ys a nS 
COM /Cmin/ Xmin.Ymin.Stx.Stry 


PRINTER IS ] 


poe 

Pei SING S2Xk 0" Reranl? valiess 3" 

D=2. (ese-3 ! Tupe diameter (m) 

Oe U ome nicht (m) 

R=,.4021 ee GU CMO mms Ont Onan 

PRIN] USING "4X," jupe Giameter mili ys) rama Cir.) cue eeres 
PRINT USING “4X.""Tube iength Ee OOD im) ha 
Sere oiNG “4k. Racius of totor = “.Z, 3p, ered (Tr) Solan 
BES» 

eo». ‘GK TG ACCEPT DEFAULT VALUES (1=Y,0=N)?".Id 

HUPUT “LIKE A BARD COPY ¢'=¥.C=N)?" ihe 

TF The=1 THEN 

PRS Meese Olt 

ELSE 

Reh Estas 

END IF 

Tr Ed={} THEN 1300 

BEEP 

THPUT “ENTER TUBE DIAMETER ¢m)".D 

BEEP 

LN SSE isis (We ENG? im) e gil 

Bee 

INPUT “ENTER RADIUS GF ROTOR (mo”".R 

SR SUNG sot OxK. a (GSO aS meters sip mo 
Serily US EME He AX) MOV sleigecer Coy = ee i) ea aaemet (ry) cna) 8) 
PRINT USING “14x. “ube rengtn (L) = DE. DD. aoscemi Cir )snvaeanerts a 
Seo NlGem hone ae WD = OD One) 
PRINT USING “14X.° "Rotor radius (R) = we 730.0" arr) Rese 1 
ce 

ORT SN Goenka Sono lice Cl wc tOSm inom mana 

panrer Sn : 

SRINT Sen Gime etna oan =e Oa me) et Neca 

YSOTIN IWIN YaNC MP tT  SiitiaKe MS exeyiintgs 

SMT STNG MA Pee Min wepesuns led tiie 

PRIN? USING "4X.°°3 Wu versus Re” 

PASI USOUNKG SC Sa IN, exeireSieSy [seven 

INPUT Go 

iF Op>3S THEN 

PRINT USING "4x, “NSELEC™ (al Tara (0) hice 

PRINT USING v4X5,- 2 VARY RPM 

COMO SNG mea xe mes | VARY Rew 

INPUT Opx 

Tf Oox=2 THEN Qn=2 

IF Gpx=3 THEN Gp=3 

END IF 

IF Up>1 THEN 

BEEP 


J:=0 

Set II EI Rm tomnt 
PRINT USING “4X. 
PRINT USING “4X, 


pooner ere iy (ake xian ee ee 
we ella Semaclioe Irae : 


2=Nakayama. 3=Nakayama & Fuziola 


106 





FOMDINMTS 
eo 
COND 
oo0o 


PRINT USING “4X,."" 4=Nooas & Morr:s. S=Stephenson. S=hoodsaMorrisze, iam 
ose vem co nr! 
INPUT Tan 


Lote 


MNrg) Like fu eLrut vizt.u=N)7°,Ukolet 
IF Ovplot=1 THEN 
pees PLC ECS F pH othals hg Wo Ue niGred 
Eom 
INPUT “SOLID=!.DASH-DASH=2 ,.DOT-DOT=3". Type 
ENS ie 
SN ee 
Nstep=100 
TF Up=2 THEN 
BEEP 
INPUT “ENTER RPM RANGE (MIN.MAX)?” .Roml.Romh 
Rom=Rpm!l 
ELSE 
BEEP 
INPUT “ENTER RPM OF ROTGR".Rpm 
END IF 
BEEP 
TNPUT “ENTER OPTION (1=0,2=Mf.3=Dt)". Io 
IF Ihc=1 THEN PRINTER IS 701 
TF Thc=0 TREN PRINTER IS 1 
PRINT 
PRINT USING “0X, === Uperating Variables *=~ ~~ 
PRINT 
Eat it NG men aes DIL ar lap lesmanhe ssn 
IF To=1 THEN 3 
BEEP 
TNPUT “ENTER CODLANT MASS FLOW RATE (ka/s)" MF 
BEEP 
OREN USING aN “Coolant Nas sen Lowa Le 
INPUT “ENTER INLET AND QUTLET TEMPS (C)".T 
9RINT USING "14X,.""Coolant inlet teme 
PRINT USING “14X.""Coolant outlet temp 
END IF 
TF To=2 THEN 
See 
INPUT “ENTER HEAT LGAD (H)".Q 
eer USING "14X."""heat load cs TONED) [DIDS PCI POM oT 
BEEP 
INPUT “ENTER INLET AND QUTLET TEMPS (C)".T 
PRINT USING “14X.""*Coolant iniet temp 
PRINT USING "14X,"""Coolant outlet temp 
END IF 
IF fo=3 THEN 
Ir J1>0 THEN 2200 


acres. (egy ape Me 


er oua OD CO ST a 
200 0D CC) 316 


noite W 
o 
o 


10 
sedan) [IN Des oe Oe OTE 
OOS (Oke OD eae GG) maeeeniC) 


Le | 


BEEP 

INPUT VENTER HEAT LOAD ¢H)", G 

PRINT USING "14X,""Heat load Se PAR Ob (CHD 
BEEP 


IF Qp<3 THEN 

INPUT “ENTER COOLANT MASS FLOW RATE (kg/s)" MF 

BREN USING “W4,° “Coolant mess Gilet) wees S 75 oS Ceeyey re atthe 
ELSE 

INPUT “ENTER MASS FLOK RATES (MIN.MAX)” .Mf1l.MFh 

MF=mfl 


Pe eeeNieR COOLANT INLST TEMP (Co”". 1: 
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= 


-- 


PRINT USING “i4X.""Coo.ant iniet temp eee UU a cea 


END iF 
hash l io. 5 
Cp=tNUpu(Ta) 
Nu=FNMNuw( Tad 
K=FNKw( Ta) 
Dr=FNPru (ia) 
Rho=FNRhow( Ta) 


Nuv=Mu/Rho 

IF Io=! THEN 

G=Mf=Cp«(To-1i) 

PRINT USING "14X,""Heat load (@Q) ame’ 40D. ace se 

BNO TF 

Ir Io=2Z THEN 

Mr =0/(Cp*(To-T1)? 

BrclieUSiNG “4X. "Coolant mass ¢low rate = °".Z.4DE,°" (kg/eo°'" MT 


IF Io=3 THEN 

Toc=Ti+Q@/(Mf Co) 

TE ABS(To-Toc)>.01 THEN 

To=(Tot+Toc)*.5 

GOTO 2220 

END IF 

IF Jy=0 THEN PRINT USING "14X.""Coolant outlet temp ee FD BID GP CE 


END IF 

IF J1=0 THEN 

PRINT 

Retin Sth Gam OG eC mp mo pemtte Sma cella tie clit D Diy) Ol am mG) etch Orient era: 


Ct USN Gamal x ee Specit NG weata Gop) See hae CK gk) Sep 
PRINT USING “i4X.""Viscosity (Cu) Se Aree) See aene ANU SOT We mae 
PRINT USING “*4X."“"Tnerma! cond (k) hearer 11) Pacem GLY AiCTy ped ae aa ace 
CN US UNG ILaXnmone GanGc Uanlimben Cra). SO ED OP ie 

PRINT USING “14X."“Density (Rho) = 000 ZN) 5D), PP Chaya SNP Bins 
PRINT USING “14X.""Kinematic vis (Nuv) S000 ADE OO ane Ble ee sii 
Beta=FNBeta(Ta) 

PRINT USING “14X.""Coef ther exp (Beta) 0 Te ADE OC EAR DOO Serer 
PRIN! 

MeotieuUSENGe: 10X.0=== Catcrlations ===) 

PRINT 

PRINT USING “10X.""8reliminary calculations:”"”" 

END IF 


Re=4=Nf/(OT=*D*Mu) 


Friction Factors for stationary reference 
IF Re<2*1.&+4 THEN 2670 

rao Ly Rem mirc) maonmonre:| 

GOTO 2530 

Borcnici hose mu LoDeeann oo 


Qpp=Q/(PI+D=L) 

Qp=Opp=PI«D 

YUm=4=MFf/(Rho*PI<D~2) 

Des=F *Rho*Vm* 2*L/(2*D) 

Umeaqa=Rpm=2*P1/60 

Ro=Vm/(Omega*D) 

Jay=Omega«D°2/Nuv 'From STEPH-Rotational Re 
Kw="NKw( 7a) 

Dtt=(To-T1)/L 
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Owonrnurwn— 
ooocoooocoonod 


69) 60 09 29 0 GI OO LO 0) 
(\) = =a + 2 2 ee ee ot 


! Nakayama Correlation for tu 


Grak=(jmega 2)-Setae((D/2) 4)<Dtr/Nuv 2 
Ra=GrePr 
Gamma=(Re © Geo Vo me CGP Tr SD mee ReaD.) 


=). 2 ee od 
TS- Ne vuecrcuantal © sc) 


! Laminar Value 


Nulam=48/11 


Turbulent Value from pees Boelter 
Nudb=.023*Re* .8=Pr- 


! Steonenson Correiation 


B@e-«'r .4/.367)*.007i*Re .88eJay .023 'Air,Turb-E0N2.23-corrected for wat 


urpulent canditions 
Nunmak=(Re>.8)*(Pr° .4)*(7Tq° (1/30) )*(1+.014/( Tg) (1/6))=.033 !'! Eqn 2.9 (Tur 


Ty=Nunak/(Re~.8*Pr .4) 


! Nakayama Correlation for laminar conditions 


J2=J~2 

Sq=2/114(1+SGR01+77/4*(1/Pr°2))) 

Cf=1-.486*((3*Sq-1)° .4/(Sq*(Sq*Pr*SOR(5)+2))>*J2/(RaxRe)’ .6 
Nunak !=48/11+.191/Sa*(3*Sq-1)7> .2*(Ra*Re)* .2*1/¢1+¢1/10*#Saq*Pr)) 


! Nakayama/Fuzioka Correlation for Radial P:pes 


Nunf=(.014/.023)=(Re/Ro°2.5)° .124*Nudb '!TURRB-EQN 2.21 
Mus=Mu 'Temp Value Nus®Mu 


' Siecer-jiate Correlation for Turbulent 


Nust=.02/7*Re’° .8=Pr* .3333=(Mu/Nus) .14 
Prod=Ra-re-Pr 


Wenoocds-Norrs Correlation tor Laminar 


Nuum=.262=Prod” .173*48/711 ' Ean 2.22-WM2 
Ji=zJay/8 


WNooas-Morris Correlation ror Radia! Yipes 
Nuym2=.015*Re® .78*J1°.25 'from wm2 
IF Js=0 THEN 


PRINT USING “14X.""Mean fiuid vel (Vm) eee (n/s00"" <Um 


PRINT USING "14X,."""Reynolds number (Re) = aya) Cae 
PRINT USING "14X.°"Friction factor (stat) = ee. 
PRINT USING “"14X,.""Heat flux (Opp) ce OO AIDE OP CRIA) > RE Pa Blea 
PRINT USING “14X,°"“RPM =e [aD maesin Cr 
END IF 
IF The=0 THEN PRINTER IS i 
PRINT 
feo OmiaaNmowrN MUSING: NOX, “oReswl ts: see: 
IF Okeplot=0 GR (Okplot=1 AND Nstep<11) THEN 
PRON USING OX. Nudb Nunak Munf Nuwm Nus Nuum2'"***"* 
PRINT USING “13X,.6(3D.DD.2X)" :Nudb.Nunak .Nunf .Nuwm.Nus,Nuwm2 
PRICN MUSING 14x, °° Ra Ro RasRe*Pr Mf (hea kgs 
PRINT USING “14X.4(Z.2DE.2X) .3D.DD.DD.DDD.Z.DD":Ra.Ro.Prod,Mf ,Nunak 1 
PRINT UISHiNGe 14x, 28" Rare Tg ily enone 
aio 1 Do ADE oso Za SWE « By, Z.DDD" :RaxRe.Tg. Ty 
N 


IF Okplot=1 THEN 
IF lan=1 THEN Y=Nudb 
IF Tan=2 THEN Y=Nunak 
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THEN Y=Nunf 

THEN Y=Nuwm 

THEN Y=Nus 

THER YeNouun? 

THEN Y=Nulam 

THEN Y=Nunak } 

IF Ge=2 AND Oox<>2 THEN X=Rpm 


-e 
nn 
tre 
u 
5) 
CON MUD Lo 


IF Dpx=2 OR Upx=3 THEN X=Prod 
Sea a tal 

Siem otG jiwG.IdaNn + lYOe. Gx) 
END IF 

TF Op=2 THEN 

Rom=Rom*10° (Cx/Nstep? 

IF Rpm>Romh THEN 3590 

GOTQ 2730 

END IF 

IF Op=3 THEN 
Mf=Mf=«10°(Cx/Nstep) 

IF Mf>oMfh THEN 3590 

GOTO 2040 

END IF 

BEEP 

INPUT “ANOTHER RUN (1=Y,0=N)?" Ir 
PRINT “PuU" 

TF Ir=! THEN 1250 

INPUT “WANT TO LABEL?(1=Y.0=N)".11 
iF Li=! THEN CALL Labei 

END 





Ee 2, ae iy) 

DIM K(8> 

HATS ot S3° 324554 ,-25 pe TS Ie = =e SOCRAE 24 Be eo RS nd =e te iad (Eee 
UM MOM Se All SiaGOnwlet wels ac 

READ Kc) 

bCisteamter3.15)/64- .2 

Sum=0 

FOR N=0 TQ & 

Sum=Sum+K(N)*C1-T) CN=1) 

NEXT N 


Br=oum/ (1<¢ tek (5) = 1-7 4K 0G eC 1-7) 72) - C1 -T)/ OK 07) #01 -T) 24K 08) dD 
Pr=cxXPc2r) 

P=22120000*Fr 

RETURN P 

SNEND 

DEF FNHFe(T) 
Hfa=2477200-2450=«(7-10) 

RETURN Hic 

FNEND 

DEF FNMuw(T) 

9=247 .38/(T+133.15) 
Mu=2.4E-5-10°A 

RETURN Mu 

FNEND 

DEF FNVvst(Tt) 

P=FNPyst(Tt) 

nt tess 

X=1500/T 

Fi=t/cl+Text ,£-4) 
F2=(1-EXP(-X))72.S5=EX P(X) /X" 5 
B= .0015<*F '- .000942=F 2- .0004882<x 
K=2*P/(461 ,52<T) 
YeCt+(1+2sh=K)° 5) /K 

RETURN V 

FNEND 

DEF FNCow(T) 

cope dee Ve aoe l= (2 72oacoe-ool (loc oollaoet oe. / d2 97 ~ 1) 
RETURN Cow=1000 

FNEND 

DEF FNRhow(T) 

Roe Comes ot (One been e 2 ooe oes oe eon) 
RETURN Ro 

-NEND 

DEF FN rut T) 
Pru=FNCow(T)=FNMuw( 7) /FNKu¢T) 
RETURN Pru 

-NEND 

DEF FNKw(T) 

KeIGiaeesenl Ie ons 

Kw=- .92247+X*(2.8395-xK<(1.8007-X*(.52577- .07344*xX))) 
RETURN Ky 

FNEND 

DEF FNTann(X) 

P=FXP(X) 

Q=EXP(-X) 

Tann=(P-Q)/¢(P+Q) 

RETURN Tanh 

-NEND 

DEF FNHf (7) 


iad 





sf=T<(4,203849-7#(5.28132E-4-7<4,5516035176£-9)) 
RETURN Hf=1000 

=NEND 

DEF FNTwerne Pd 

yu-110 

= 

Ta=(Tutil<.5 

Po=FNPvst (ja) 

T& ABS CP -Po)/P)>.001 THEN 

IF Pce<P THEN Tl=Ta 

IF °¢>? THEN TueTa 

GOTa i650 

=ND TF 

RETURN Ta 

FNEND 

DEF FNBeta(t) 

Rop=FNRhow( T+. 1) 

Rom=FNRhow(T-.1) 

Beta=-2/(Roept+Ram) =(Rore-Rom)/.2 
RETURN Beta 

FNEND 

DEF FNAlpha(T) 
Alpha="NKu(T)/¢rFNRhow(T)*FNCow(T)) 
RETURN Alpha 

“NEND 

DEF FNSigma€T> 

7t=547.3-T 
A=.001<Tt*2*¢.1160936807/(1+.83*Tt)) 
B=.001121404688-5.75280518E-6-Tt 
Cares ee boa ore aS ioe Sites 
Srama=A+5+C 

RETURN Sigma s 

FNEND 





+7 Seane 


metric Yarreplec <== 
Tune diameter (D) = 4,753E-03 (nm) 
Ree on seh es = Boe 
7A = er 
Roter radius (R) = 0.402 (m) 
Computed values are: 
-*> Jperating Variables <= 
iD viGcmny inset n ere eee a: 
Heals = $9.00 Ci) 
Coclant :nlet temp = 45,00 (C) 
Coolant outlet temp CS e/a ie 
Coolant mass flow rate = 1.6329E-02 (kg/s) 
SUicm ROO im tNe sme Vv Hella to cm sttme Mo om CDmone!, 
‘Specific neat (Cp) Se eee ee Fees ce) eK) 
Viscosity (Mu) = 5.8661£-04 (N.s/m°2) 
Thermal cond (k) = 0.6380 (W/m.¥) 
Prandt! number (Pr) = 3.884 
Density (Rho) = 990.0 (kg/m~3) 
Kinematic vis (Nuv) =O Acca Om CMa cis) 
Coet ther exp (eeta) = 4,1266E-04 ¢1/K) 
=e Calculations **<* 
Preliminary calcuiations: 
Reynolds numper (Ke) 7. 445,£+03 
EBTpciee TO men GIG OGM Gite cttcy) mee ines eg il ene th 
Heat flux (Gop) = 4,060EF+03 (W/m?) 
Mean fluid vel (Ym) = 9.26£-0! (m/s) 
REM = 3600.00 
Resicse 
Nudb Nunak Munt Nuwm Mus Nuwim2 
Gorse Some 417.89 31.56 44.84 102.34 
Ra RO Rarka<Pr i4¢ Nak] 
(Eos GROG mn Oa Cael oat an Ee ESO 44,90 
Rake SG Tv 
S.S0&+07 7.506£-05 0.026 
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APPENDIX C 
R 


THERMOSYPHCN ANALYSIS (COMPUTER PROGRAM AND RESULTS) 


The following proaqram was written to give the results of 


Reomeco-c-tst.om for the heat transfer for the rotating, 
closed-lcop thermosyfhon discussed in Chapter Iil. Ti was 
Metweeh =meHP-Basic for tha HP-9826 computer. Sic o-ograme 
Be= tne Graphics and the chermophysical proverties are not 
shown. A sample run is included immediately follewing the 
tise To. 


-Pages 114-117 Program Listing 
rage 118 Sample Calculation 





2a me 


=a MAME: 


CREMICY! 


CQM /Cain/ X 
PRINTER IS 3 
Soar 

PRINT USING 
D=4.7635-3 
e=no250 
R=.4021 
KFe=Q 

ce eee 
PRINT USING 
PRINT USING 
PRINT USING 
PRINT USING 
BEEP 


ce 


“SpUrarVy ¢0.. one 


min, ¥min.57x,Stv 


Sho” Waganulte wellness 
' Tube ciameter (m) 
' “ube iength Ca) 
Pacis Om nO LO mm CmD 
anicuLonmeacton Ton 


"4X.°"Tupe diameter 
"4X.""Tupe Lengtn 
"AX." "Ckt Length (Lt) 
"4X .°" Radius of rotor 


bends 


Soe aDe m 


tous .OD . D be 
ee DD 
soon We 63D mest 


Meu ok TO ACCEPT SEPRUET VALUES (1=Y.0=N)?" Id 


BEEP 


Mee 2=ike A APRD COPY «t=¥,0=N0?" The 


Pe ine=1 THE 
Tr id=] THEN 
ED 


iNPUT “ENTER 
BEEP 
INPUT "ENTER 
BaP 
PAPUT “ENTER 
PRINT USING 


PRINT USING ' 


PRINT USING 
PRT USING 
PRINT USING 


PRINT USING ° 


PRINT 

PRINT USING 
PRINTER IS 7 
BEEP 

PRINT USING 
PRINT USING 
PRINT USING 
INSUT Op 

Tr Go>i THEN 
Baar 

Jys=9 
ReneS. | 


NEPRINTER=1S6 701 
1280 


TUBE DIAMETER (m>".D 
TUBE SLENGIH <a)". L 
RADIUS OF ROTOR (m)".R 


"10X,°"=*+ Geometric Variables ~** 
Saar 


BOX ee OCMC iam elte na Dn) 
ate een Cm eT Git mG le 

ae -1) ramus (Ocak CH at re G emt) 
eee , ark ait 

HAY PNG eee ieelel Ae: (RD) 


WiUKP ee econputed «alies anes 


poe oe Ci eIP TION 


BY ci amcce” | Single Dont saan 


cee ae ene Ulan eee 


Pegs sne 1G PLOT (1=7.0=N)?7" .Okplot 


BEEP 


2DE.°" 


sees ,DD OD mete 
tous JOD OD a 
Paola seer D) ' 


le 


Ae 


INPUT “Nudb=1.Nunak 1=2,.Nuum=3.Nus=4.Nul=5"..Tan 


IF Okplot=! 
CALL Plot(Jy; 
NP ieee Ke 
IF Il=1 THEN 
END IF 

END IF 
Nstep=11 

iF Op=2 THEN 
BEEP 


INPUT “ENTER REM RANGE (MIN. MAX)" Romi .Romn 


THEN 

piso tt o Daan iteo) 

A LABEL7(1=Y.0=N)".11 
CALL Lapel 


yoveeee Dd 


Cayes oi 
Cy eee Se 
Cron tk 


Gm 
Cm aie 


Gn 


Greet: 


sj 
. 


-* 


ai 


D 


+ 





1650 
1500 
1670 
1630 
1699 
*700 
1710 
1720 
1730 


750 


i780 
cele 
13800 
B10 


LNFUL CEMTOR APM fr tri tIR™ Dan 


IF The=1 THEN PRINTER IS 7n1 
if Inc=0 THEN PRINTER IS 3 

PRINT 

PRINT USING “!0X.°"=** Qperatina Variables *** °°" 
PRINT 

Bie ot GCmmnUA sent van iasles jared = 

BEEP 

TNR oR BAD CW) .G 

PRINT USING “14X,""Heat load = meee) ea) eae 
cop 

aneUieeenMen CUED SHDS TEMP (CC). ie 
PRINT USING "14X.""Cold side temo 
Dt=5 

Th=Tct+Dt 

a=(To+th)*.5 

Dine De /e 

PRINT 

Co=FNCouw(Ta? 

Mu=FNMuw( 7a) 

K=FNKw( Ta) 

Pr=FNPrwtTa) 

Rho=FNRhow( Ta) 

Nuv=Mu/Rho 

Beta=FNBeta(ia) 

Omega=hom=2*P1/60 

Th=Te+Dt 

Ta=ic+Dt/2 

PRINTER IS 7 

PRINT “TA = “Ta 

Rhoh=FNRnow( Th) 

PRINT "RHOH="',Rnon 

Rhoc=FNRhow( Tc) 

PRINT “RHOC=".Rhoc 

Rhom=FNRhouw( Ta) 

Drho=Rhoc-Rhon 

IF Drho<O THEN GOTO 3:40 

PRN  acdro=..0 rho 

Mu=F NMuw (Ta) 

Co=FNCow( Ta) 

Beta=FNBeta(Ta) 

Pr=FNPru(Ta) 

Nuv=Mu/Rhom 

Prati UinGe s=canMin, Co 

Mf=0/(Cpe*Dt) 

Re=4/Pl«MfF/(Mu*D) 
Gr=R*(Omega*2)*Beta*((D/2)° 4)*Dit/Nuv 2 
Ra=Gre*Pr 

Prod=Ra*Re*Pr 
Nuwm=.262*Prod~.173*48/11 
Dtwo=0/(K*Nuwm*P IL) 

PRINT “DTHP =",Dtup 

Dtw=Dtwp 

Grr=R«Umega”2<Beta*D 3*Dtw/Nuv 2 
Nu=Pr/GrrxLt/R=D/L*Rhom/Rhon* (Kf F+(256/Re*Lt/D))*Ke 2 
Dtwc=G/(PI«L<k«Nu) 


aoee 4D .DD as 


116 


Gib) emcee a 2) 


CGD ae aie 





TF ABS (Dtwe-Dtud> 70 THEN 
TF Rpm>3000 THEN 
Nesh~+s C 


VUTU 1 oUY 


END IF 

IF Rpm>!000 THEN 

Dt=Dt+} 

GOTG 1900 

END IF 

IF Rom>S00 THEN 

Die= Dit .> 

GOTO +S00 

END IF 

TF Rom>=300 THEN 

Dt=Dtr+4 

GOTO :3900 

END IF 

END IF 

Ym=4*MF/(Rhom=P I<) 2) 

Ro=Vm/(Omega+D) 

Jay=Omega+D 2/Nuv 'From STEPH-Rotationa! Re 
Det=(Th-Te)/L 

Gr=R*(Gmega 2)-Beta=((D/2) 4)*Dtt/Nuw "2 
Ra=Gre-Pr 


Laminar Value 
Nulam=48/1}3 


Turbulent Value from Dittus-Boelter 
Nudb=.023*Re .8*Pr’ .4 


Stephenson Correiation 
Nus=(Pr° .4/.367)-.007!=Re .o8-say ’.023 !Air. Turb-cON2.23-corrected for wat 


ie Correlatiagn tor laminar cond:it:ons 

=e? 

Sq=2/15<(1+SORCi+77/4-(1/Pr -2))) 

Cf=1-.486<((3*Saq-1)° .4/(Sa*(SqePr=SQR(5)+2)))*J2/(RasRke) .6 
Nunak 1=48/11=.191/Sa<(3#Sq-1)° .2=(RasRe) 281/01 +(4/10-+Saq~Pr)) 
Mus=Mu 'Temp Value Mus®Mu 

Prod=Ra*Re#Pr 


Hoods-Morris Correlation for Laminar 
Nuwm=.262*Prod* .173*48/11 ! Ean 2.22-WMS 
Jl=Jay/8 


Woods-Morris Correlation for Radial Pipes 

Nuwm2=.015<Re* .78*J1°.25 'from wm2 

IF the=0 THEN 

PRINTER IS 1 

elLSE 

PRINTER IS 701 

=<ND IF 

IF Jy=0 THEN 

IF Jy=0 TREN PRINT USING “10X.""Results: ”“"" 

ClNimUSINGmlOX Qe ritircmorcoer tle smevalttatecma tes DDD Dl mn GG) cimelnis rs 


PRAT USMIG “14X.""Seec:fic heat (Co? nae Et) Carman GUI Ka) ime SO 
PRINT USING “14X.""Viscosity (Mu) MN Z. ADE CNS sm coe Mu 


“tt 


ay 





-_— 


WWW OIG) IDI) 0) 
WONIMWN Swi 
ooooco00OGO 


ee wk oe et we et 


PRINT USING “14X.""Temperature (hotside) 


sO UUs ree 
PRONT USING “T4X.° "RPM 


PQEINT HSING “tax, herma!l cand (i) Saas eer ol fi} enue NE TT] Pep Ca) ae neeat 
SOR GUS G cay. totam Smo e Gacas eo 5 Syokes | Be 

PRINT USING weay i “Density (rho? =e error he Kec Amo! Rho 
~~ EM Ty tr RIO Mae cc, se Ste sie Cees: es 
Snevi UorNomir1k, “Goer therm axp’neca) SO eae a AO IR meta 
PRINT USING “14X,""Reynolds number (Re) = we 5 apes Re 

PRINT USING "14X."""Mean fluid ve! (Vm) SE il Ula Po aA 
PRINT USING “14x, ass flow rate  Eaeien 7) ey tate 


© DD0DO.. <Rhom 


END IF 

PRINT 

fF Gkoios=0 OR (Gkpiot=" ANd pater sto THEN 

PRINT USING “10X."" 9 Nudb = Nunak! = Nuwm Nus 0" 
PRINT USING °73X.6(3D.DD, 2X)" :Nudb .Nunak ! .Nuwm.Nus 

PRINT USING “14X,""' Ra Ro RatkexPr Mf reese 
PRINT USING “14X,4(Z.2DE.2X).2D.DD":Ra.Ro.Prod.Mf 

END IF 


TF Qkplot=1 THEN 

IF Tan=1 THEN Y=Nudb 

IF Ian=2 THEN Y=Nunak | 
IF Tan=3 THEN Y=Nuuwm 

IF Ilan=4 THEN Y=Nus 

IF Tan=5 THEN Y=(48/11) 


Uae 

GA mia a tiGuunexe Ye Dite. le) 

PRINTER IS i 

END Ie 

IF Gpe=2 THEN 

Dt=5 

Rom=Rom<10°.1° 

Tf Rpm>Rpmh THEN 3130 

GOTG i890 

END Ir 

BEEP 

INPUT “SNOQTHER RUN (4=Y.0=N)?" [ir 
PRINT “Pu 

IF Ir=! THEN 3360 

INPUT “WANT TO LABEL?¢(1=Y.Q=N)".11 
IF ITl=! THEN CALL VTapel 

END 





metric 


Tupe diameter 


Yarcsioae 


(D>? 


Mou 


as 
ft 
bata 


Hoe We Cae 


Tube length (L? = Seam) 
Ckt length (Lip 2.45 (m) 
L/D = Wea 

Rotor radius (R) = 9.402 (m) 

Computed values ara: 

See Ons anim cdma oc Smee © 

otihevcantaples are: 
saeac load = HCG) 

Coia siae teme = Seer mn oy) 

Results: 

Fluig proverties evaluated at 32.75 (C) are: 
Specific neat (Cp) = 4310.1E+00¢J/bae.k) 
Viscosity (Mu) ae erage ce Cal) ence) 
Tnermal cond (k) = 0.5404 (W/m.k) 
Prandtl number (Pr) 2.302E+00 
Density (Rno) eee mmm race) 
<inematic vis (Nuv) PSS hia eres) 
Coef therm exp( Beta) Sens Aiea) clea) 
Reynolds numper (Re) = |,978E+02 
Mean fluid vel (Vm) =e 4 ee ema) 
Mass flow rate = 3.27£-04 (ka/s) 
Temperature (hotside) = 8C.50(C) 

RPM ~ = 3600. 
Nudp Nunal | Nuwm Nus 
eee SoG: 34.44 aS 
Ra Ri Ra-Re<Pr fe 
Sree See On 2 ek 02 ee pie et) merase! 


Lis 





ADPENDIX 9 
Gree2D TYO-PHASE THERZOSYPHON BNALYSIS (COMPUTER PROGRAM AND 
RESULTS) 


The following pregram was written to give the results of 
Cer 


the EelatlOons f£O2 the «sonic, boiling, and condensing 
diwi=s as discussed in Chapter IV. [a ews werteon in 
HP-Basic for the AP-9826 computer. Sub-programs for the 


graphics and the thermophysical properties are not shown. A 
ample run is included immediately following the listing. 


"9 


ages 120-124 Program Listing 
-Pages 125-127 Sample Calculation 


120 





“a 


ce ~ 
e 
1930 


‘oayt 


19350 
19690 
0970 
7980 
1030 
10 


2 Sy 
eo an 
eed 
1120 
1749 
msc 
50 
ee 0 
119 
ieee) 
1200 
oye 
Be 
.a5 
co 
“240 
Teo) 
1250 
he 


i440 
1450 
1460 
1470 
1430 
1490 
1500 
1519 
520 
io 0 
1540 
$350 
+566 
ef) 


"580 


ia. 


ft 





ee =Ea7 

Sian hee a3: 
Bees oe Mar - 72, 394 
Se etrecos, S000 920; S250 FNS: 6700 ,caRe. 
eed ne An. Sn ot <.oty 
Sores mens 
1 Tee 6 toe aaa 
SBae 
Seka Gee Nc eran = valiess vite 
D=4./53=-3 ! Tupe diameter (m?) 
ease Ae } jing Lemmy (Svar ¢ a) 
ene Pes 
Spas 
22302" ' “eqgius a> sator (mn) 

eee 3 et . arate aie 
SES ENG 94 (DR diameter = nner eae) =a Cm) cat 
Bee elgeG) “2, upe lengrntEvand= “".90.0.°" ae 
POT NS ENG cs emer Tee ~anathncs and)= oeee ie ae eeee ( beseee ars 
See SMG “IX. 2Saqius of rotor = "",7.2D.°"" cay nenS 
Dea) 
a hen ae 
EMPYS “OK FR ACCEST SESAUILT VALUES (1=¥ J] 2". Te 
Bes 

ec ON eM ap Cle =r] OC ae 
Be Paes > Fac 
ee et) 
ELGS 
PRINTER IS 
EMD IF 

Pao ee ob 0 
Boz? 
RI eS SE DEANERER (nm)? .D 
BEE? 
INPUT “ENTER TUBE LENGTH(Evap) (m)". Le 
pase 
TMI “Sls Sey Tle Wena G(Geray Cm gle 


3E=> 
TMPUT 
PRINT 
PRINT 


aera 
SING 
“SSING 


RADIUS OF 
ite 


s+ 
ane ;4X, 


ROTOR 
Geometric 
“’Tupe diameter 


Cm)" 
Var:abies 
«D) = eee 


ere eee 


Ae we Cm)" 


PRINT USING re ““Tupe jiength(Evap) (Le) MeL OU: 
Perit USiNG “!4¥.°"Tube lenatn(Cond) (Lc) = “".Z.30.°" 
PVG Geno XwmeROcor radius sR) eae ee Den 


P| Ti 
PRINTER 
BEEP 

PRINT 
PRINT 
OQ ENT 
INPUT 


a 
& 
Le 


seueee 


USING 
USING 
USING 
Sp. 


Lee Sots ti ie Eel 
ae waey 
aeaeth Ker acer 


DP oes 


ier SUcsiee ane ea 
DT MCD atten meee 


IF Jp>i THEM 
pees 
POINTER IS 1 


Pasty USeiG “4X 
GUT 144G 
Se 


= 
cS 


SONOT IMPLEMENTED YET" 


Best 
5+) 
deat. 
25s 

BAe) 


aan Ge eC ie) =") amen Oy NOre 


aie 


OTD) 
Cine fh 
Coe 2! 


Cm pee oR 


= 


ce 





ee - 
SYS SAMOS NM MT Jas TFSi n 


MA EOIFs - OMIM EOIRs 


ery 


ee 


a yey 
Spd ede 
“20 

SS 


a ~~ eee wt 
2 =e) 
o3 


Z° 30 


eh Goede Gass ap Cra te ad act 


~aoOen 


oS 


Patel 
: a 
oO 


T= Sy == = 
we ee emo seas, to Ee 
deco 
ea Naa to 
Seen Syl ol, = Re Simo ia 
eNO I= 
END 1° 
Nstan=100 
Bea 
HIP Sagan TP pein Ch=S) 2s 
[© ioe=2 THEN 
SaaS 
tego ee = SS VER, fe 
CMO Lt 
Soe eo aayires 2S 784 
TF The=0 THEN PRINTER IS ° 
Se 
ie eee eee cca) mane a 
sete WASH IN Gammaane ZX fame loci Can 
aeee 


INPUT “ENTER RPM" .Rpn 
SRINT JSING “Tax ,"" REM 


Dimes 
SHPOT “ENTER CEOMD. WALL TEMPER 
Ree Mme mmm Ch Xeni tye! 


Ace °t<) 2/4 
cain =is ciate ©. 
Soa eomiel sco) 
Mul= "NMuwt (Tsa 
Mave" ‘Miu tsat 
K="=NKw(Tsat) 
Pis=epiPsiuis tsa) 
Rhoi=-NRhow(Tsat) 
Hi q="NRFalTsat) 


us 


ania. 
) 


Beta=FNBetal saz? 
Stama="‘IS:ama( sat) 
Rhove’ /ENVvstlisat) 
Nuv=4iuw/Rhow 

RSI =n eer 

Cv=Co/Rsn 
Ree = Solio cum GU kaj ike) 
Umega="pm<*2<PT/60 
Accel=Qmeaa 2<% 

LF a) THEM 

PRIN” 


Pe Sn Git) AP tomo nGooonmtme smo 7alllWwate oillait 


ARTS USSG OF Bi, MOSES RF he AYES 


PRINT USTHG “14X,.""Rario of 
PEO NGmmA nt EE SGOs ll cy 


PLORE 7% 


(Co) 
Seecific 
(Mu) 


PRINT USING “14X Thermal cond ¢r) 


ITTY WIS OMG Be, Masini inidinexsie 
(nf q) 


PRINT USEHE “14X.“"Entaalpy 


Kn GenGur An savensitey ORhol) 
PRINT USING "14X%.""Deneity (Rhov) 


pues USSG “l4X.""Coet tner exp (Be 


PRINT USING "14X.""Omega 


PRINT BS LG eam hi Tal mC SOln 


EMD. sh 


ICl = om 
m ~? 


122 


(Pr) 


pl 
oon uw 


tous O07 
. 


a see CNY 
Dee Dae i 
“7 ADE. ee ch 
Sea Se we ae) ser Sek 
FADE" Heo 

i COR am SY sini 
pO Ra igy BDO Be 
» os NOE PP CAR ee Blain ee, 
Art [Doe ene Od nee 28 TN ey 
GAD aS SES cl 


Pea cls 


POD OO hu 


—U00RO0. "(Caras 


maine ae ‘Co 


s/m eo" sMul 





74 care: 
— : a . 
yan 
—--v¥ 


221M Av=(M/C_t-(Rhov-8hai)))-C (Pho L/(RHOv-Rhol))+At 
aaah 


2230! SGNTC SALCULSTION 

2240 TF To=1 THEN 

Sao Mee TSE; 

2260 PRINT 

mee SP SIG “29K. "Sonie Limit Calculations:""" 

228 BRU SENG NOX ye 4 eal Gites amt ona amen 

27290 END IF 

700 Ome Hhowereqe(4/ (Ute (Rhou-Shol))-Rhol /(Rhow- Rhol)«At)<SGR(CRs ee ote laine) 
2345) MP=iim/ rts 


Mano 9 be See Pp en Bret 


Boe | fa2 5 

2230 Y=Qm 

2320) f= Bm<=9 FEN GETO 2400 
2350 fF Okpiot=? TREN 


POCO en Peioc CU eG yoCR A ny Oe.6 x? 
Zee Se 

2730 Sal) icra; Gee ees Dino Nees eGo em 
eee etl Dia 


Bay) Seo ce >= "90 Te Garg 3720 
efi) fF Okelot=i THEN 

ete Persees 1) 02 

2a) giles 

2440 Per=ePet-} 

B85 Ira aye! 

Beco) Sal) ke 

2am “Gear 220C 

2430) SEs! £P 

24539! 


Za0n! Sth aoe ‘ENT CALCULATICNS 


IT 
poe te lo=2 THEN 


Bee. 
eeou ap ie) WEN 


2540 PRINT 

SoU ee ews GMOS ants ning mimi t Graciat vons: 
2550 SPRINT USING “70X58 °% Fil: lotal Hage ss: 

eof) ID LF 

See0' MAPLE IS CORRELATION J<=.25 

2590 IF Te=‘ THEN 

2500 IF 4v=) THEN GOTO 2 
2610 Jv=.25/Rhov .S<(Ome 
2620 Om=Khoavedvedeedtg 
2630 END IF 

2540! 


Q 
“2*R<D*(Rhoi-Rhov))” .5 


ao 
ega 


2630' JASTER CORRELATION 

2550 IF [e=2 THEN 

2570 ir Pets) THEN Pst=l 
2680 T1=€4/¢1-(PeF/1N0)))-1 
25390! Ti=t/(Pef/100))-! 

2700 T2=71<¢(Rhol/Rhov)~ (2/3) 
2710 TF T1<0 THEN 2990 


2729 eeu! /Muw) .25<1! 1.754 (Rhol /RRov? > €(2/3)<1.75-1) 
2730 Xx=1/¢(4+T2) 

PM X= 

2750! <x2=(hul/Muv)  .25<2hov/Rho! 

2750 IF G=0 THEN G=: 

ees ee 46-0 /MuL) 

27380 Rev=aX%y+G-D/Muv 


oi) =v=.0/79/F%euv 1.25 


2 


23 





) 
oD 
oO 


t 
u 
¢ 


oww win 
omnawc 


) 


? 
ad 


WOhuprol se 
@iorinroe 


= 


Fain) 
= 2 
(ve) 
a) 
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